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Summary 

A detailed exposition on a refined nonlinear shell theory that is suitable for nonlinear limit- 
point buckling analyses of practical laminated-composite aerospace structures is presented. This 
shell theory includes the classical nonlinear shell theory attributed to Leonard, Sanders, Koiter, 
and Budiansky as an explicit proper subset that is obtained directly by neglecting all quantities 
associated with higher-order effects such as transverse-shearing deformation. This approach is 
used in order to leverage the exisiting experience base and to make the theory attractive to 
industry. In addition, the formalism of general tensors is avoided in order to expose the details 
needed to fully understand and use the theory in a process leading ultimately to vehicle 
certification. 

The shell theory presented is constructed around a set of strain-displacement relations that are 
based on "small" strains and "moderate" rotations. No shell-thinness approximations involving 
the ratio of the maximum thickness to the minimum radius of curvature are used and, as a result, 
the strain-displacement relations are exact within the presumptions of "small" strains and 
"moderate" rotations. To faciliate physical insight, these strain-displacement relations are 
presented in terms of the linear reference- surface strains, rotations, and changes in curvature and 
twist that appear in the classical "best" first-approximation linear shell theory attributed to 
Sanders, Koiter, and Budiansky. The effects of transverse-shearing deformations are included in 
the strain-displacement relations and kinematic equations by using analyst-defined functions to 
describe the through-the-thickness distributions of transverse-shearing strains. This approach 
yields a wide range of flexibility to the analyst when confronted with new structural 
configurations and the need to analyze both global and local response phenomena, and it enables 
a building-block approach to analysis. The theory also uses the three-dimensional elasticity form 
of the internal virtual work to obtain the symmetrical effective stress resulants that appear in 
classical nonlinear shell theory attributed to Leonard, Sanders, Koiter, and Budiansky. The 
principle of virtual work, including "live" pressure effects, and the surface divergence theorem 
are used to obtain the nonlinear equilibrium equations and boundary conditions. 

A key element of the shell theory presented herein is the treatment of the constitutive 
equations, which include thermal effects. The constitutive equations for laminated-composite 
shells are derived without using any shell-thinness approximations, and simplified forms and 
special cases are discussed that include the use of layerwise zigzag kinematics. In addition, the 
effects of shell-thinness approximations on the constitutive equations are presented. It is 
noteworthy that none of the shell-thinness approximations appear outside of the constitutive 
equations, which are inherently approximate. Lastly, the effects of "small" initial geometric 
imperfections are introduced in a relatively simple manner, and a resume' of the fundamental 
equations are given in an appendix. Overall, a hierarchy of shell theories that are amenable to the 
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prediction of global and local responses and to the development of generic design technology are 
obtained in a detailed and unified manner. 


Major Symbols 


The primary symbols used in the present study are given as follows. 




A 

A,, A, 


unit-magnitude base vector fields of the shell reference surface 
shown in figure 1 

area of shell reference surface (see figure 1), in 2 

metric coefficients of the shell reference surface 
(see equation (1) 


A-115 A 12, ^ 16? A 2 2,A 26 , 


shell stiffnesses defined by equation (101) 


An, A12, Ai 6 , A 2 2, A 2 6, A 66 
A 44 , A 45 , A55 

Bn, Bj 2 , B 16 , B 22 , B 2 6, B 66 

Cp C 2 , C 3 

C„ 

[C u ] 

ds 

Dll? D12, Dig, D 2 2, D 2 6, D 66 


[d 0 ], [d,], [d 2 ] 


shell membrane stiffnesses (see equation (B31)), lb/in. 
shell transverse-shearing stiffnesses (see equation (B32)), lb/in. 

shell coupling stiffnesses (see equation (B31)), lb 

tracers used to indentify various shell theories 
(see equation (53)) 

transformed shear stiffnesses appearing in equation (85b), psi 

constitutive matrices defined by equations (87) 

differential arc length define by equation (1), in. 

shell bending and twisting stiffnesses (see equation 
(B31)), in.-lb 

matrices defined by equations (26), (27), and (55) 


[do], [d\], [d;] 


matrices defined by equations (11 6) - (11 7) 


000 
^11? ^22? ^12 


linear deformation parameters defined by equations (8) 


3 



Fj(jy , F,(! 3 ) analyst-defined functions that specifiy the through-the-thickness 

distributions of the transverse-shearing strains (see 
equations (3)), in. 

?n^ 12 ,? 2l ,? 22 A?)A? 2 } work-conjugate stress resultants defined by equations (20c) 

and (20d), lb 

G(|,) function defining the through-the-thickness temperature 

variation (see equation (90)) 

gn, gt g£ shell thermal coefficients defined by equation (109) 


h 


shell-wall thickness, in. 




maximum shell thickness divided by the minimum principal 
radius of curvature 


h„, h ]2 , h 2 


shell thermal coefficients defined by equation (108) 


K’ k 45> k 


55 


[k 0 ], [k,], [k 2 ], [kj, [k 12 ], [k 22 ] 


M n , M 12 , M 21 , M 2 




M,(y, M 12 (y 


M 21 (^j), M 2 (|j) 


transverse-shear correction factors appearing in equation (B32) 

matrices defined by equations (28) and (29) 

bending stress resultants defined by equations (13), in-lb/in. 

work-conjugate bending stress resultants defined by equation 
(20b), in-lb/in. 

applied loads on edge = constant (see figure 2), in-lb/in. 
applied loads on edge '§ 2 = constant (see figure 2), in-lb/in. 


unit-magnitude base vector field perpendicular to the shell 
reference surface, as depicted in figure 1 


N 

N N N N 

Ml’ x 1 12’ - l ^21 9 x '22 




N,0y 


n 2 (U 


unit-magnitude vector field perpendicular to the shell reference- 
surface boundary curve and n , as depicted in figure 1 

membrane stress resultants defined by equations (13), lb/in. 

work-conjugate membrane stress resultants defined by equation 
(20a), lb/in. 

applied load on edge = constant (see figure 2), lb/in. 
applied load on edge t; 2 = constant (see figure 2), lb/in. 
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V v V v P 3 

p\. Pi p\ 


q,, q 2 , q 3 , {q} 


D D D 

Q 1 ? ^2 ? ^3 


q 3 


Ql? ^2? ^3 


Q,(y 




Ql3’ Q 23 


■^ 13 ’ ^ 23 ’ {^} 

Q.3, Q 23 


Qy 


Qu, QS, Qt QL Q*\ Qs6 


effective tractions defined by equations (59), psi 

effective tractions defined by equations ( 122 ), psi 

applied surface tractions (see equations (32)), psi 

dead-load part of applied surface tractions (see equations 
(32)), psi 

live-load part of applied surface tractions (see equations 
(32)), psi 

tractions associated with interactions between live pressure and 
initial geometric imprefections, defined by equation ( 122 d), psi 

applied load on edge = constant (see figure 2 ), lb/in. 

applied load on edge = constant (see figure 2 ), lb/in. 

transverse-shear stress resultants defined by equations (13), 
lb/in. 

work-conjugate transverse- shear stress resultants defined by 
equation ( 20 e), lb/in. 

stress resultants defined by equations (59), lb/in. 

transformed, reduced (plane stress) stiffnesses of classical 

laminated- shell and laminated-plate theories 
(see equation (85a)), psi 

shell stiffnesses defined by equation (103) 


R p R principal radii of curvature of the shell reference surface along 

the and < 5 , coordinate directions, respectively, in. 


R ' k R jk R ' k R jk R jk R 

1 1 5 xv 12’ 16’ rv 22’ rv 265 1X 


jk 

66 


s.cy 


S 2 (^) 


[S 0 ], [S,], [SJ, [S 3 ], [S 4 ], [S 5 ] 


shell stiffnesses defined by equation ( 102 ) 

applied load on edge = constant (see figure 2 ), lb/in. 

applied load on edge = constant (see figure 2 ), lb/in. 

matrices defined by equations (17) 
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Up u 2 , u 3 
Up U 2 , U 3 


w‘fe, W 


■Jk .jk ijk 

S^-44 , -ZV45 , -z*-55 


displacements of material points comprizing the shell reference 
surface (see equations (3)), in. 

displacements of shell material points (see equations 
(3)), in. 

known, measured or assumed, distribution of reference-surface 
initial geometric imperfections measured along a vector normal 
to the reference surface at a given point, in. 


shell stiffnesses defined by equation (104) 


Y' Jk y' Jk Y' Jk 

A 44 9 A 45 9 A 55 


shell stiffnesses defined by equation (105) 


Z P Z 2 


Z 


7 uk 7' jk 7' jk 

a -'44 9 *-‘45 9 -^55 


quantities defined as i + — and l + — , respectively, and 

R, R 2 

used in equations (89) and (99) 

quantity defined as Zj + z 2 + ^(z 2 - z, ) and used in 
equations (89) and (99) 

shell stiffnesses defined by equation (106) 


a 


0 


transformed coefficients of thermal expansion appeaing in 
equation (85a), °F 1 


Yl3> Y23, (y"} 

Up 


transverse- shearing strains evaluated at the shell reference 
surface (see equation (16c)) 

transverse shear function defined by equation (5g) 


Se°p be° 22 , 6e° 2 variations of the linear deformation parameters defined 

by equations (23) 

6e n , 6 e 22 , (V/ i2 , (V/ i3 , 5y 23 , 6s 33 virtual strains appearing in equation (14b) 

„ O o o / °\ 

o£ ip o£ 22 , oy 12 , y o£ j virtual membrane strains defined by equation (22a) and (53) 

SYi 3 > by 23, {by } virtual transverse-shearing strains appearing in equation (22c) 

6cpp 6qp 2 , bcp virtual rotations of the shell reference surface about the 

and '§,-axcs, respectively, defined by equations (23), radians 


6 



fiXn’SxL. sxL {sx°} 


vector of virtual bending strains defined by equation (22b), in 1 


m 

ml 


m x 


m 

ml 

ml, ml 


external virtual work per unit area of shell reference surface 
defined by equation (32a), lb/in. 

external virtual work per unit length of the applied tractions 

acting on the boundary curve dA that encloses the region A 
(see figure 1 and equation (33a)), lb 

external virtual work per unit length of shell reference surface 
boundary defined by equations (33), lb 

internal virtual work per unit area of shell reference surface 
defined by equations (19), lb/in. 

internal virtual work per unit volume of shell defined by 
equations (14), psi 

internal virtual work per unit length of shell reference surface 
boundary defined by equation (47), lb 

internal virtual work per unit length of shell reference surface 
boundary defined by equations (48), lb 


SUj, 6u„ 6 u 3 , {8u} 


virtual displacements of the shell reference surface about the 
| r , and § 3 -axes, respectively, in. (see equation (27a)) 


£ 11? e 22> Y 12; Y 13? Y23> £ 33 


shell strains defined by equations (5) 


£ in 


Yn? {e"} 


reference-surface normal and shearing strains defined by 
equations (6) and (51) 


K n , K 22 , 2 k 12 changes in reference surface curvature and torsion defined by 

equations (7), in 1 

p u , p 22 radii of geodesic curvature of the shell reference surface 

coordinate curves and respectively, in. 
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o,i, o 2 21 a u , a 13 , o 2} , o 33 


T 



«Pp <P 2 > <P 


shell stresses, psi 

parameter used to identify second-order terms in equations (100) 

curvilinear coordinates of the shell, as depicted in figure 1 

linear rotation parameters for the shell reference surface defined 
by equations (4), radians 


° ° O ° / ° \ 

2 X,»[Xj linear deformation parameters defined by equations (8) 

and (16b), in.' 1 


0,0 


function describeing the pointwise change in temperature from 
a uniform reference state (see equations (90)), °F 


{0 k }, {©J 

BA 


thermal quantities defined by equations (88) and (92), 
respectively 

curve bounding area of shell reference surface (see figure 1), in. 


Introduction 

Classical plate and shell theories have played an important role in the design of high- 
performance aerospace structures for many years. As a result, familiarity with these theories is 
generally widespread throughout the aerospace industry, and a great amount of resources has gone 
into validating their use in design. Perhaps the best-known classical shell theory is the one 
attributed to A. E. H. Love 1-2 that was re-derived by Reissner 3 in 1941. This shell theory has 
become known worldwide as the Love-Kirchhoff classical thin-shell theory. This particular shell 
theory, as presented by Reissner, has some deficiencies that were later addressed by Sanders, 
Budiansky, and Koiter 46 to obtain what is generally deemed as the "best" first-approximation 
classical thin-shell theory. This shell theory was later extended to include the effects of geometric 
nonlinearities by Leonard, 7 Sanders, 8 Koiter, 9 10 and Budiansky. 11 For the most part, these theories 
are focused on shells made of isotropic materials. 

As the need for improved structural performance has increased, new materials and design 
concepts have emerged that require refined plate and shell theories in order to predict adequately 
the structural behavior. For example, a sandwich plate with fiber-reinforced face plates and a 
relatively flexible core, either of which may have embedded electromechanical actuation layers, 
is a structure that typically requires a refined theory. Similarly, efforts made over the last 20 to 30 
years to reduce structural weight or to enable active shape control have resulted in thin-walled, 
relatively flexible designs that require nonlinear theories to predict accurately responses such as 
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buckling and flutter. 

Many refined plate and shell theories have been developed over the past 50 or so years that 
are classified as equivalent single-layer, layer-wise, zigzag, and variational asymptotic theories. 
Detailed historical accounts of these theories are beyond the scope of the present study and can 
be found in references 12-80. Each of these theories has its own merits and range of validity 
associated with a given class of problems, and the choice of which theory to use depends generally 
on the nature of the response characteristics of interest. For the most part, these theories have not 
yet found wide acceptance in standard industry design practices because of the extensive 
experience base with classical theories, the relatively limited amount of validation studies, and the 
increased complexity that designers usually try to avoid. In general, validation studies associated 
with stuctures made from exotic state-of-the-art materials are very expensive if experiments are 
involved. Moreover, there are usually many more structural parameters that must be examined in 
order to understand the design space, compared to the number of parameters that characterize the 
behaviors of the more commonplace metallic structures. 

The present study is concerned with the development of refined shell theories that include the 
classical shell theories as well-defined, explicit proper subsets. Herein, the term "explicit proper 
subset" means that the equations of a particular classical shell theory appear directly when all 
quantities associated with higher-order effects, such as transverse shearing deformations, are 
neglected. In contrast, the terms "implicit subset" and "contained implicitly" are used to indicate 
cases where the equations of a particular classical shell theory can be recovered by using a 
transformation of the fundamental unknown response functions. This interest in refined shell 
theories that include the classical shell theories as well-defined, explicit subsets is motivated by 
the need for design-technology and certification technology development that takes full 
advantage of the existing experience base. For example, legacy codes used by industry that have 
undergone extensive, expensive experimental validation over many years can be enhanced to 
address issues associated with new materials and design concepts with a high degree of 
confidence. Moreover, this approach appears to avoid undesirable computational ill-conditioning 
effects. 61 Fikewise, experience and insight gained in the development and use of nondimensional 
parameters 81165 to characterize the very broad response spectrum of laminated-composite plates 
and shells can be retained and extended with the high degree of confidence needed to design and 
certify aerospace vehicles. Furthermore, the development and use of nondimensional parameters 
have a high potential to impact the development of scaling technologies that can be used to design 
sub-scale experiments for validation of new analysis methods and for flight certification of 
aerospace vehicles (e. g., see reference 147). 

Of the many refined theories for plates and shells discussed in references 12-80, several are 
particularly relevant to the present study. 166 ' 212 In an early 1958 paper by Ambartsumian, 166 a 
general equivalent single-layer, linear theory of anisotropic shells was derived that presumes 
parabolic through-the-thickness distributions for the transverse-shearing stresses. Subsequent 
integration of the corresponding strain-displacement relations is shown to yield expressions for 
the displacement fields that include those of the classical theory of shells explicitly as a proper 
subset. Six equilibrium equations are also used that involve the asymmetrical shearing and 
twisting stress resultants that are obtained by integrating the shearing stresses across the shell 
thickness. A similar derivation was presented later by Ambartsumian 167 for shallow shells in 1960. 
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Likewise, Tomashevski 168 used a similar approach to derive the corresponding equations for 
buckling of orthotropic cylinders in 1966. 

In 1969, Cappelli et.al. 169 presented equations for orthotropic shells of revolution that are 
based on Sanders’ 8 linear shell theory and that include the effects of transverse-shear 
deformations. In those equations, the two rotations of a material line element that is perpendicular 
to the shell reference surface are used as fundamental unknowns and, as a result, the 
corresponding equations of classical shell theory do not appear explicitly as a proper subset. In 
contrast, Bhimaraddi 170 presented linear equations for vibration analysis of isotropic circular 
cylindrical shells in 1984 that include parabolic through-the-thickness distributions for the 
transverse- shearing stresses and contain Flugge’s equations 171 as an explicit proper subset. It is 
noteworthy to recall that Flugge’s equations retains terms of second order in the ratio of the 
maximum thickness to the minimum radius of curvature that is used in the shell-thinness 
approximations. 

Also in 1984, Reddy 172 presented a linear first-order transverse-shear-deformation theory for 
doubly curved, laminated-composite shells that extends Sanders’ original work 4 by including the 
two rotations of a material line element that is perpendicular to the shell reference surface that are 
used as fundamental unknowns and by introducing constitutive equations that relate the 
transverse- shear stress resultant to the transverse shearing strains. Similarly, in 1985, Reddy and 
Liu 173 extended Reddy’s previous shear-deformation theory for doubly curved, laminated- 
composite shells by including parabolic through-the-thickness distributions for the transverse- 
shearing stresses. Like Cappelli et.al., 169 the equations in references 172 and 173 do not contain 
the equations of Sanders’ shell theory as an explicit proper subset. 

Soldatos 174 178 presented a refined shear-deformation theory for isotropic and laminated- 
composite non-circular cylindrical shells during 1986-1992. This particular theory includes 
parabolic through-the-thickness distributions for the transverse-shearing stresses and contains the 
equations of Love-Kirchoff classical shell theory as an explicit proper subset. Additionally, only 
five independent unknown functions are present in the kinematic equations, like first-order 
transverse-shear deformation theories. In 1989, Bhimaraddi et al. 179 presented a derivation for a 
shear-deformable shell finite element that is based on kinematics that include parabolic through- 
the-thickness distributions for the transverse-shearing strains in addition to the kinematics based 
on the hypothesis originally used by Love. 1 ' 2 Likewise, in 1992, Touratier 180 presented a 
generalization of the theories discussed herein so far that combines parabolic through-the- 
thickness distributions for the transverse-shearing strains and the classical Love-Kirchhoff linear 
shell theory. Specifically, following his earlier work on plates (see reference 181), Touratier 
appended the Love-Kirchhoff kinematics for shells undergoing axisymmetric deformations with 
a transverse-shear deformation term that uses a somewhat arbitrary function of the through-the- 
thickness coordinate to define the distributions of the transverse-shearing stresses. The 
arbitrariness of this function is limited by the requirement that the corresponding transverse 
shearing stresses satisfy the traction boundary conditions on the bounding surfaces of the shell. 
This process yields general functional representations for the transverse shearing stresses, much 
like that of Ambartsumian, 166 167 182 that are specified by the analyst a priori. Moreover, by 
specifying the appropriate shear-deformation functions, the first-order and refined theories that 
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have uniform and parabolic shear-stress distributions, respectively, discussed herein previously 
are obtained as special cases. Touratier 181,183 also presented results based on using sinusoidal 
through-the-thickness distributions for the transverse-shearing strains that are similar to those 
used earlier by Stein and Jegley. 184,185 A similar formulation for shallow shells was given by 
Sklepus 186 in 1996, which includes thermal effects. Additionally, in 1992, Soldatos 187 presented a 
refined shear-deformation theory for circular cylindrical shells that is similar to general 
formulation of Touratier 180 but utilizes only four unknown functions in the kinematic equations 
and also accounts for transverse normal strains. Later, in 1999, Lam et.al. 188 determined the 
vibration modes of thick laminated-composite cylindrical shells by using a refined theory that 
includes parabolic through-the-thickness distributions for the transverse-shearing strains in 
addition to the kinematics of classical Love-Kirchhoff shell theory. In 2001, Fares & Youssif 189 
derived an improved first-order shear-deformation nonlinear shell theory, with the Sanders-type 
kinematics used by Reddy, 172 that uses a mixed variational principle to obtain stresses that are 
continuous across the shell thickness. A similar theory was also derived by Zenkour and Fares 190 
in 2001 for laminated cylindrical shells. 

Recently, Mantari et.al. 191 presented a linear theory for doubly curved shallow shells, made of 
laminated-composite materials, that is similar in form to the derivation given by Reddy and Liu, 173 
but uses the form of the kinematics used by Touratier 180 in 1992 and discussed previously herein. 
In contrast to Touratier’ s work, the theory given by Mantari et.al. contains the linear equations of 
the Donnell-Mushtari- Vlasov 192 shell theory as an explicit proper subset. In addition, Mantari et 
al. use a special form of the functions used to specify the through-the-thickness distributions of 
the transverse-shearing strains that contains a "tuning" parameter. This parameter is selected to 
maximize the transverse flexibility of a given laminate construction. A similar derivation, but 
with an emphasis on a different form of the functions used to specify the through-the-thickness 
distributions of the transverse-shearing strains was presented by Mantari et al. 193,194 in 2012. Very 
recently, Viola et.al. 195 presented a general high-order, linear, equivalent single-layer shear- 
deformation theory for shells that contains many of the theories described herein previously as 
special cases. 

Several shell theories have been derived over the past 25 to 30 years that utilize layerwise 
kinematics to enhance an equivalent single-layer theory without introducing additional unknown 
independent functions that lead to boundary- value problems of higher order. In 1991, Librescu & 
Schmidt 196 presented a general theory of shells that appends the kinematics of first-order shear- 
deformation shell theory with layerwise functions that are selected to yield displacement and 
stress continuity at layer interfaces. However, the traction boundary conditions at the top and 
bottom shell surfaces are not satisfied. In addition, the theory includes the effects of relatively 
small-magnitude geometric nonlinearities. In contrast, to a standard first-order shear-deformation 
shell theory, this theory of Librescu & Schmidt has a twelveth-order system of equations 
governing the response. 

Later, in 1993 and 1995, Soldatos and Timarci 197,198 presented, and applied, a general 
formulation for cylindrical laminated-composite shells, similar to that of Touratier, 180 that includes 
five independent unknown functions in the kinematic equations and a discussion about 
incorporating layerwise zigzag kinematics into the functions used to specify the through-the- 
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thickness distribution of transverse shearing stresses. Similarly, in 1993 and 1995, Jing and 
Tzeng 199-200 presented, and applied, a refined theory for laminated-composite shells that is based 
on the kinematics of first-order transverse-shear-deformation shell theory appended with zigzag 
layer-displacement functions, and on assumed independent transverse-shearing stress fields that 
satisfy traction continuity at the layer interfaces. A mixed variational approach is used to obtain 
the compatibility equations for transverse-shearing deformations in addition to the equilibrium 
equations and boundary conditions. Although, the theory captures the layerwise deformations and 
stresses, it has only seven unknown functions in the kinematic equations, regardless of the number 
of layers. In addition, the theory includes the exact form of the shell curvature terms appearing 
in the strain-displacement relations of elasticity theory and in the usual, general definitions of the 
stress resultants for shells. Likewise, in 1993, a general theory for doubly curved laminated- 
composite shallow shells was presented by Beakou and Touratier 201 that incorporates zigzag layer- 
displacement functions and that has only five unknown functions in the kinematic equations, 
regardless of the number of layers. Similar work was presented by Ossadzow, Muller, Touratier, 
and Faye 202 203 in 1995. Moreover, Shaw and Gosling 204 extended the theory of Beakou and 
Touratier 201 in 2011 to include non-shallow, deep shells. 

In 1994, Fie 205 presented a general linear theory of laminated-composite shells that has the 
three unknown reference-surface displacement fields of classical Love-Kirchhoff shell theory and 
two additional ones that are selected to satisfy continuity of displacements and stresses at layer 
interfaces and the traction boundary conditions at the two bounding surfaces of a shell. Similarly, 
Shu 206 presented a linear theory for laminated-composite shallow shells in 1996 that also satisfies 
continuity of displacements and stresses at layer interfaces and the traction boundary conditions 
at the two bounding surfaces of a shell. Shu’s theory, however, contains the equations of the 
Donnell-Mushtari- Vlasov 192 shell theory as an explicit proper subset. In 1997, Shu 207 extended this 
theory to include nonshallow shells, with the classical Love-Kirchhoff shell theory as an explicit 
proper subset. 

Cho, Kim, and Kim 208 presented a refined theory for laminated-composite shells in 1996 that 
is based on the kinematics of first-order transverse-shear-deformation shell theory appended with 
zigzag layer-displacement functions. In contrast to the theory of Jing and Tzeng, 199 this theory uses 
a displacement formulation to enforce traction continuity at the layer interfaces. In 1999, Soldatos 
& Shu 209 presented a stress analysis method for doubly curved laminated shells that is based on 
their earlier work (e.g., see references 197, 198, 206, and 208), which includes five unknown 
independent functions in the kinematic equations and two general functions that are used to 
specify the through-the-thickness distribution of transverse shearing stresses. The two functions 
that are used to specify the through-the-thickness distribution of transverse shearing stresses are 
determined by applying the two equilibrium equations of elasticity theory that relate the 
transverse-shearing stresses to the stresses acting in the tangent plane to obtain a system of 
ordinary differential equations for the two unknown functions. This approach yields solutions for 
the two functions in each shell layer. The constants of integration are determined by enforcing 
continuity of displacements and stresses at layer interfaces and the traction boundary conditions 
at the two bounding surfaces of a shell. 

The present study is also concerned primarily with the development of refined nonlinear shell 
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theories. Several previous works relevant to the present study are given by references 210-220 and 
72. Specifically, in 1987, Librescu 210 presented a general theory for geometrically perfect, elastic, 
anisotropic, multilayer shells of general shape, using the formalism of general tensors, that utilizes 
a mixed variational approach to obtain the equations governing the shell response. These 
equations include continuity conditions for stresses and displacements at the layer interfaces. In 
his theory, the displacement fields are expanded in power series with respect to the through-the- 
thickness shell coordinate, and then substituted into the three-dimensional, nonlinear Green- 
Lagrange strains of elasticity theory. This step yields nonlinear strain-displacement relations with 
no restrictions placed on the size of the displacement gradients, and it contains the classical "small 
finite deflection" theory of Koiter 9 and the classical "small strain-moderate rotation" theory of 
Sanders 8 as special cases. A wide range of refined geometrically nonlinear shell theories can be 
obtained from Librescu ’s general formulation, each of which is based on the number of terms 
retained in the power series expansions. In 1988, Librescu and Schmidt 211 presented a similar 
derivation for a general theory of shells, again using the formalism of general tensors, that uses 
Hamilton’s variational principle to derive the equations of motion and boundary conditions. This 
particular derivation did not yield continuity conditions for stresses and displacements at the layer 
interfaces, and focused on geometric nonlinearities associated with "small" strains and 
"moderate" rotations. Likewise, Schmidt and Reddy 212 presented a general first-order shear- 
deformation theory for elastic, geometrically perfect anisotropic shells in 1988, following an 
approach similar to Librescu and Schmidt, for "small" strains and "moderate" rotations that 
includes uniform through-the-thickness normal strain. Their derivation also uses the formalism of 
general tensors. Another similar presentation and an assessment of the theory was given by 
Palmerio, Reddy, and Schmidt in 1990. 213,214 In contrast to the theory of Librescu and Schmidt, 
the theory of Schmidt and Reddy utilizes a simpler set of strain-displacement relations which 
neglects nonlinear rotations about the vector field normal to the reference surface. Moreover, the 
dynamic version of the principle of virtual displacements is used to obtain the corresponding 
equations of motion and boundary conditions. Furthermore, the classical "small" strain and 
"moderate" rotation theories given by Leonard, 7 Sanders, 8 and Koiter 9 10 are contained in the 
Schmidt-Reddy theory implicitly; that is, they can be obtain by using a change of independent 
variables in the equations governing the response. 

In 1991, Carrera 215 presented a first-order shear-deformation theory for buckling and 
vibration of doubly curved laminated-composite shells that includes the Fliigge-Lur 'e-Byrne 
equations 216 as an explicit proper subset. This particular set of classical equations for doubly 
curved shells also retains terms of second order in the ratio of the maximum thickness to the 
minimum radius of curvature that is used in the shell-thinness approximations. The geometric 
nonlinearities used by Carrera are identical to those of the Donnell-Mushtari- Vlasov shell theory, 
as presented by Sanders. 8 In 1992 and 1993, Simitses and Anastasiadis 217 ' 218 presented a refined 
nonlinear theory for moderately thick, laminated-composite, circular cylindrical shells that 
includes geometric nonlinearity associated with "small" strains and "moderate" rotations, like 
that given by Sanders, 8 and initial geometric imperfections. The theory is based on cubic through- 
the-thickness axial and circumferential displacements and constant through-the-thickness normal 
displacements, and neglects rotations about the vector field normal to the reference surface. 
Moreover, the classical theory of Sanders 8 is contained as an implicit subset. Also in 1992, 
Soldatos 219 presented a refined nonlinear theory for geometrically perfect laminated-composite 


13 


cylindrical shells with a general non-circular cross-sectional profile. Soldatos’ theory presumes 
parabolic through-the-thickness distributions for the transverse-shearing stresses and contains the 
equations of Love-Kirchoff classical shell theory as an explicit proper subset. Additionally, only 
five unknown independent functions are present in the kinematic equations, like first-order 
transverse-shear deformation theories. The geometric nonlinearity correspond to the "small" 
strains and "moderate" rotations of Sanders, 8 with rotations about the normal vector field 
neglected. The equations of motion and boundary conditions are obtained by using Hamilton’s 
variational principle. 

Several years later, in 2008, Takano 220 presented a nonlinear theory for geometrically perfect, 
anisotropic, circular cylindrical shells. His theory uses the full geometric nonlinearity possessed 
by the three-dimensional Green-Lagrange strains of elasticity theory, and the equilibrium 
equations and boundary conditions are obtained by applying the principle of virtual work. 
Moreover, the theory is formulated as a first-order shear deformation theory and, when linearized, 
includes Fltigge’s equations 171 as an explicit proper subset. Takano’s work is similar to that 
presented in 1991 by Carrera 215 for doubly curved shells, but includes a higher degree of 
nonlinearity in the strain-displacement relations. In 2009, Pirrera and Weaver 72 presented a 
nonlinear, first-order shear-deformation theory for geometrically perfect anisotropic shells that 
uses the rotations of material line elements normal to the reference surface as fundamental 
unknowns. In their theory, the full geometric nonlinearity possessed by the three-dimensional 
Green-Lagrange strains of elasticity theory is used and expressed in terms of linear strains and 
rotations, and their products. Additionally, the equations of motion used in their theory are based 
on momentum balance of a differential shell element, as opposed to being determined from a 
variational principle. Moreover, the equations of motion are linear, which appears to be 
inconsistent with a geometrically nonlinear theory. 

The literature review given previously herein reveals a need for a detailed exposition on a 
refined nonlinear shell theory that is suitable for nonlinear limit-point buckling analyses of 
practical aerospace structures made of laminated-composites that utilize advanced structural 
design concepts. A major goal of the present study is to supply this exposition. Another goal is to 
focus on a shell theory that includes the classical nonlinear shell theories as explicit proper subsets 
in order to leverage the exisiting experience base and to make the theory attractive to industry. To 
accomplish these goals, the formalism of general tensors is avoided in order to expose the details 
needed to fully understand and use the theory in a process leading ultimately to vehicle 
certification. In addition, the analysis is simplified greatly by focusing on the many practical cases 
that can be addressed by using principal-curvature coordinates. The key to accomplishing these 
goals is the form of the strain-displacement relations. 

The strain-displacement relations used in the present study are a subset of those derived in 
reference 221, that are useful for nonlinear limit-point buckling analyses. These strain- 
displacement relations are based on "small" strains and "moderate" rotations, and are presented 
first, along with a description of the shell geometry and kinematics. Moreover, the strain- 
displacement relations are presented in terms of the linear reference-surface strains, rotations, and 
changes in curvature and twist that appear in the classical "best" first-approximation linear shell 
theory attributed to Sanders 4 , Koiter 5 , and Budiansky. 6 The effects of transverse-shearing 
deformations are included in the strain-displacement relations and kinematic equations by using 
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the approach of Touratier 180 in which the through-the-thickness distributions of transverse- 
shearing strains are represented by analyst-defined functions. Additionally, no shell-thinness 
approximations involving the ratio of the maximum thickness to the minimum radius of curvature 
are used and, as a result, the strain-displacement relations are exact within the presumptions of 
"small" strains and "moderate" rotations. This approach yields a wide range of flexibility to the 
analyst when confronted with new structural configurations and the need to analyze both global 
and local response phenomena. Next, the usual asymmetrical shell stress resultants that are 
obtained by integrating the stresses across the shell thickness are defined, and the three- 
dimensional elasticity form of the internal virtual work is given and used to obtain the 
corresponding symmetrical effective stress resulants that appear in classical nonlinear shell theory 
attributed to Leonard, 7 Sanders, 8 Koiter, 910 and Budiansky. 11 Afterward, the principle of virtual 
work, including "live" pressure effects, and the surface divergence theorem are used to obtain the 
nonlinear equilibrium equations and boundary conditions. Then, the thermoelastic constitutive 
equations for laminated-composite shells are derived without using any shell-thinness 
approximations. Simplified forms and special cases of the constitutive equations are also 
discussed that include the use of layerwise zigzag kinematics. In addition, the effects of shell- 
thinness approximations on the constitutive equations are presented. It is noteworthy to mention 
that none of the shell-thinness approximations discussed in the present study appear outside of the 
constitutive equations, which are inherently approximate due to the fact that their specification 
requires experimentally determined quantities that are often not known precisely. Lastly, the 
effects of "small" initial geometric imperfections are introduced in a relatively simple manner, 
and a resume' of the fundamental equations are given in an appendix. Overall, a hierarchy of shell 
theories are obtained in a detailed and unified manner that are amenable to the prediction of global 
and local responses and to the development of generic design technology. 


Geometry and Coordinate Systems 

The equations governing the nonlinear deformations of doubly curved shells are presented 
subsequently in terms of the orthogonal, principal-curvature, curvilinear coordinates (§,,^ 2 ,§ 3 ) that 
are depicted in figure 1 for a generic shell reference surface A. Associated with each point p of 
the reference surface, with coordinates (^,,^,,0), are three perpendicular, unit-magnitude vector 
fields a,, a„ and n. The vectors a, and a 2 are tangent to the and ^-coordinate curves, 
respectively, and reside in the tangent plane at the point p. The vector n is tangent to the in- 
coordinate curve at point p and perpendicular to the tangent plane. The metric coefficients of the 
reference surface, also known as coefficients of the first fundamental form, are denoted by the 
functions A^,,^) and A 2 (|,,| 2 ) that appear in the equation 

(ds) 2 = (A,d^) 2 + (A 2 d^ 2 ) 2 (1) 


where ds is the differential arc length between two infinitesimally neighboring points of the 
surface, p and q. This class of parametric coordinates permits substantial simplification of the 
shell equations and has many practical applications. 
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Principal-curvature coordinates form an orthogonal coordinate mesh and are identified by 
examining how the vectors a,, a 2 , and n change as the coordinate curves are traversed by an 
infinitesimal amount. In particular, at every point q that is infinitesimally close to point p there is 
another set of vectors a,, a,, and n with similar attributes; that is, the vectors a, and a, are 
orthogonal and tangent to the and ^-coordinate curves at q, respectively, and reside in the 
tangent plane at the point q. Likewise, vector n is tangent to the ^-coordinate curve at point q 
and perpendicular to the tangent plane at point q. Next, consider the finite portion of the tangent 
plane at point p shown in figure 1. Because of the identical properties of the vectors a p a,, and 
n at every point of the surface, an identical, corresponding planar region exists at point q. 
Therefore, the vectors a,, a 2 , and n at point q can be obtained by moving the vectors a,, a,, 
and n at point p to point q. In addition, the plane region at point p moves into coincidence with 
the corresponding plane region at point q as the surface is traversed from point p to point q. 
During this process, the plane region at point p undergoes roll, pitch, and yaw (rotation about the 
normal line to the surface) motions. The roll and pitch motions are caused by surface twist 
(torsion) and curvature, respectively. The yaw motion is associated with the geodesic curvature 
of the surface curve traversed in going from point p to q. When a principal-curvature coordinate 
curve is traversed in going from point p to q, the planar region at point p undergoes only pitch 
and yaw motions as it moves into coincidence with the corresponding region at point q. Rolling 
motion associated with local surface torsion does not occur. This attribute simplifies greatly the 
mathematics involved in deriving a shell theory. 

In the shell-theory equations presented herein, the functions R,(^„^ 2 ) and R, ('§,,'§,) denote 
the principal radii of curvature of the shell reference surface along the and coordinate 
directions, respectively. Similarly, the functions p n (^,J§ 2 ) and p 22 (| p | 2 ) denote the radii of 
geodesic curvature of the shell reference surface coordinate curves and '% 2 , respectively. 
Discussions of these quantities are found in the books by Weatherburn 222 , Eisenhart, 223 Struik, 224 
and Kreyszig. 225 These functions are related to the metric coefficients by the equations 
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Kinematics and Strain-Displacement Relations 

The kinematics and strain-displacement relations presented in this section are special cases 
of those given in reference 221. These equations were derived on the presumption of "small" 
strains and "moderate" rotations. Moreover, no shell-thinness approximations were used in their 
derivation. In the subsequent presentations, the term "tangential" refers to quantities associated 
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with the tangent plane at a given point of the shell reference surface. In contrast, the term "normal 
refers to quantities perpendicular to the tangent plane at that given point. 


The tangential and normal displacement fields of a material point (§ 1 , ? 2 , § 3 ) of a shell are 
expressed in orthogonal principal-curvature coordinates as 


17,(1,, 1 2I E 

h) = u, + l 3 [(p, - cpcp 2 ] + F,(^ 3 )y° 3 

(3a) 

u 2 (i„i 2 ,i 

3) = u 2 + l 3 [cp 2 + qxp,] + F 2 (^ 3 )y° 3 

(3b) 

u 3 (i 

1. § 3 ) =U 3 - ^ 3 («Pl + T 2 ) 

(3c) 


where U p U 2 , and U 3 are the displacement-field components in the and ^-coordinate 

directions, respectively. The functions u,(§„ § 2 ) and u ,(§i, t 2 ) are the corresponding 


tangential displacements of the reference-surface material point (§,, i 2 . 0), and u 3 (l 1 , i 2 ) is 

the normal displacement of the material point (?i, i 2 , 0). In addition, the functions cp,(?i, i 2 ), 

tp 2 (§i» § 2 ), and <p(!i. § 2 ) are linear rotation parameters that are given in terms of the reference- 
surface displacements by 
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The functions Fj(i 3 )Y 13 (ii, 1 2 ) and F 2 (H 3 )y 23 (?i, ? 2 ) define the transverse-shearing strains. In 

particular, F 1 (| 3 ) and F 2 (| 3 ) are analyst-defined functions that specifiy the through-the- 

thickness distributions of the transverse-shearing strains, and are selected to satisfy the traction- 
free boundary conditions on the transverse-shear stresses at the bounding surfaces of the shell 

given by the coordinates ( In addition, FT?,) and F 2 (| 3 ) are selected to satisfy the 

conditions 11,(1,, 1 2 , 0 ) = u,(^, i 2 ) and U 2 (l„ 1 2 , 0 ) = u 2 (l„ 1 2 ). Thus, from equations (3) it 
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follows that F,(0) = F 2 (0) = 0. 


The nonlinear strain-displacement relations obtained from reference 221 are given as 
follows. The normal strains are 
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and the shearing strains are 
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From equations ( 5 e) and ( 5 f), it follows that Y^l^ 1 ’ ^ 2 ’ *-*) = F/(0)Yi3 and 

Y 23 (li>^2> o) = F 2 '(o)y 23 . Thus, it is convenient to scale F,(| 3 ) and F 2 (| 3 ) to give 

F,'(0) = F 2 '(0) = 1. For this scaling, Y 13 and Y 23 are the transverse-shearing strains at the shell 
reference surface. 


00 o 

In equations ( 5 ), £„> £ 22 > and Y 12 are the reference-surface normal and shearing strains, 
which are given in terms of the linear strain and rotation parameters e° lt e 22 , e° 2 , cp,, cp 2 , and cp by 
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Likewise, the changes in reference surface curvature and torsion K ,,, K 22 , and 2 k 12 caused by 
deformation are given by 
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where Xu> X22’ ar >d X12 are linear strain parameters associated with bending and twisting of the 

shell reference surface. The linear strain parameters are given in terms of the reference-surface 
displacements and linear rotations by 
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(8b) 
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In these equations, e° 1; e 22 , and 2e° 2 are recognized as the linear reference-surface strains given 

o ° 1 _ ° 

by Sanders in reference 4. Likewise, Xu’ X 22 ’ and 2 X <2 are the linear bending-strain measures 
given by Sanders. 


To arrive at the particular form of the nonlinear strain-displacement relations used in the 
present study, equations (5a), (5b), and (5d) are first re-arranged to get 
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has been used to obtain this particular form of these equations. As pointed out by Koiter, 6,910 terms 
involving a reference-surface strain divided by a principal radius of curvature are extremely small 
and can be added or neglected without significantly altering the fidelity of the strain-displacement 
relations. Thus, it follows that 
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and equations (9) reduce to 
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Equations (12) and equations (5e) and (5f) constitute the nonzero nonlinear strain displacement 
relations of the present study. The membrane reference- surface strains defined by equations (6) 
are identical to those used by Budiansky 11 and Koiter, 9, 10 and contain those used by Sanders 8 as a 
special case. Likewise, as mentioned before, the linear bending strain measures defined by 
equations (8d)-(8f) are identical to those used in the "best" first-approximation linear shell theory 
of Sanders, Budiansky, and Koiter. In contrast to the linear bending strain measures used in 
classical Love-Kirchhoff shell theory, those given by equations (8d)-(8f) vanish for rigid-body 
displacements (see reference 4). 


Stress Resultants and Virtual Work 


In the classical theories of shells, two-dimensional stress-resultant functions are used to 
represent the actual force per unit length produced by the internal stresses acting on the normal 
sections, or faces, of a shell given by constant values of the reference-surface coordinates and 
§ 2 . On an edge given by t;, = constant, the stress resultants are defined as 
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In these definitions, the middle surface of the shell is used as the reference surface, for 
convenience. Likewise, on an edge given by = constant, the stress resultants are defined as 
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where o u , o 22 , o l2 , o l3 , and o 23 are stresses and where, in general, h = h('= , , | 2 ) is the shell thickness 
at the point (| p r g 2 ) of the shell reference surface. Equations (13) show that the stress resultants 
are not symmetric; that is, N 12 ^ N,, and M p ^ M,,, even though the stresses a 12 = a,,. As a result 
of this asymmetry, constitutive equations that are based on equations (13) are typically more 
complicated that the corresponding equations for flat plates. 

To obtain symmetric stress resultant definitions that yield a simple form for the constitutive 
equations, the internal virtual work is used. The internal virtual work of a three-dimensional solid 
is given, in matrix form, in terms of the curvilinear coordinates used herein by 



where the superscript T in equation (14b) denotes matrix transposition and where A denotes the 
reference-surface area of the shell. The functions 6 e u , 6 s 22 , &y 12 , 5y 13 , 6y 23 , and 6e 33 in equation 
(14b) are the virtual strains that are obtained by taking the first variation of the corresponding shell 
strains. To obtain the form needed for the present shell theory, it is convenient to express the shell 
strains given by equations (12), (4e), and (4f) in matrix form as 
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Taking the first variation of equations (15) gives 
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Next, substituting equations (18) into equations (14) and performing the through-the-thickness 
integration yields the internal virtual work as 
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are defined as work-conjugate stress resultants. The relationship of these quantities with those 
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defined by equations (13) are 



which are identical to the effective stress results first defined by Sanders. 4 The virtual reference- 
surface strains are given by 
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where the variations of the linear deformation parameters are given by 
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' a 2 a^ 2 a, p n p 22 


S<Pl(ll.^2) = 


6 u, 1 d 6 u, 


r, a, a?, 


Scp 2 (^,, ? 2 ) = 


6 u 2 1 36u, 


R 2 A 2 d% 2 


) _ 1 

l d5u 2 

l 36u, 

5u, 

5u 2 

1 - 2 

[a, a|, 

a 2 d§ 2 

' Pn 

P22 J 


From these expressions, it follows that 


te°n = 


qh e , 2 + cp 
R. Pn 


6 u, 


1 0 

1 + e 1 


c l+e 11& j . 0 \ 1 d6u, 

Su 2 + — — Su 3 + (1 + e„)— — r^- 

10 v 'a, 


R. 


1 o , 'I 1 dbu 2 
( e i 2 + ( p)“7" 


a, 6 %, <Pi a 1 


1 d 6 u. 


6 e 22 = * + 622 6 u, + 


<Jh _ e l2 - cp 

R, P22 


c 1 + e 22 5 / o \ 1 d&Uj 

ou, h ou, + )e 12 - cp J- 


(l + e 22 ) 


R2 

1 dSu 2 1 d 6 u 3 
A 2 d % 2 ^'A, d§ 2 


a, di§ 2 


8Yk = 


5 P 2 1 +e 22 e ° 2 + (p 

R, 


Pn P22 

+ 


6u, + 


qq _ e °2 - <P _ 1 +e° 

R, P 1 1 P 20 


6uo + 


e i2 - <p e i2 + q) 

Ri R2 


6u, 


/■ o \ 1 <36u, , o \ 1 d6u, , o \ 1 

( e - " ‘p) a;^it + ( + ( 1 + c -)t 


d5u 2 


d6u, 


a 2 v a, d|, 

1 d6u, 1 dbu 3 


( e; = +tp )iif-^ ; i ; -' p 'A,a Sa 


and 


(23c) 

(23d) 

(23e) 

(23f) 

(24a) 

(24b) 

(24c) 


28 


(25a) 


* o 1 d | 1 ) 
OY,, = — — OU, 

A, a^, I R, 


+ 


Q O ^ 1 

6)C 22 = 7. I" 


P22 Ri A 2 3 | 2 lR 


1 

J 36U; 

1 6u, 1 

1 

a5u 3 

1 

d 

( 1 

dSu 

R, A, <3^ 

Pn R 2 Pn 

a 2 

a § 2 

A, 

a|. 

[a. 

a?, 

1 ) 

8u + ^ 

l a5u 2 1 

1 

d5u, 

1 

a 

1 

a6u 

R’ j 

- R 2 

A 2 ^§2 P22 

A, 

aii 

a 2 

ai 2 

A 2 

3^2 


82x°* = 


1 d 

( 1 ] 

, 1 

( 1 | l\ 

5u, + 

1 d 

1 1 

1 

( 1 1 1 y 

A 2 dh, 2 

^ R ‘ ; 

2p„ 

(R. +R :j 

A, 3?, 

l R J 

2p 22 

R, + R:j 


5m 


( 3 _ 1) 

1 1 a8u > . 11 

f 3 _ 1 ) 

1 1 d5u 2 1 | 

f 1 a8u 3 ^ 

Ir. r 2 J 

A 2 d% 2 2l 

Ir 2 r.J 

Q. 

1 

jxn 

<~o 

< 

(A, d%, j 


+ 


1 

1 a6u 3 ' 

1 d 

1 a6u 3 ' 

1 a 

1 a8u 3 '| 

P 22 

( A 2 d'% 2 i 

A 2 a§ 2 

(A, a§, J 

A.ai, 

(A, a§ 2 j 


Equation (22a) is now expressed in matrix form as 


{ 5 e°} = [d 0 ]{ 8 u) + [d,] “|y( 8u } + [ d 2] XT^|;{ 6u } 


where 


[d.] = 




fH 




( 6u ) = 


c 

r 




W 



cpi | e° 2 + cp 


1 + e 

0 

11 

1 + e n 

R| Pn 


Pn 


R, 

1 + e 22 

Cp 2 

0 

e 12 

-cp 

1 0 

1 + e 22 

P 22 

r 2 

P 

22 

R 2 

Cp 2 | 1 ^22 J ^ 12 Cp 

qpi e° 2 - cp 

1 0 

1 + e„ 

e° 2 -<p ! e° 2 + cp 

R 1 P 1 1 P 22 

R 2 

Pn 

P 22 

Ri R 2 


1 + e n 

e° 2 + cp -cp, 


[d,] = 

0 

0 

0 



e° 2 - cp 

1 + e 22 — cp 2 



(25b) 


(25c) 


(26) 


(27a) 


(27b) 


(27c) 
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[d 2 ] = 


0 0 

o 1 o 

e 12 - qp 1 + e 2 


i ° 

1 + e , 


e 12 + cp 


0 

-qp2 

- <Pi 


Likewise, equation (22b) is now expressed in matrix form as 


where 


{Sx'i = [k.]( 8 u) + [k,] ^|-{Su} + [k ;] 


-[k„L 


1 9 

A.ai, 


{ 5 u} 


[k 1 — 1 a 

L Ki2 J a 


A 2 d§ 2 


a, a§, 


{ 6 u} -[k 22 U 


A, 3§ 1 ^A 2 a§ 2 '- f 
A 2 dJ; 2 \ A,d| 


( 8u > 

/ 

\ 8u } 


[ko] 


1 


A, dli 

1 


R, 


p 22 Ri 


1 dll 


A 2 ^^2 \ R 1 


+ 


1 


1 1 

2pn 1 Ri R 2 


PnR 2 


1 3 


An d^ 7 


V , 


1 d 

( 1 ] 

1 

( 1 

A, d|. 


2p 22 

( R , 


+ ■ 


[ki] = 


1 

0 

0 


0 

0 


ii 

2 R 2 


i 


o 

l 

p 22 

l 

p^ 


(27d) 


(28) 

0 

0 (29a) 


(29b) 


30 


0 


0 


1 


[kj = 


o 


l 

r: 


2 Ir, r ,} 


o 


pn 

0 

1 

P22 


(29c) 



00 1 
000 
000 


(29(1) 


[k I2 ] 


000 
000 
00 1 


(29e) 


[ k 2 2 ] 


000 
00 1 
000 


(29f) 


Equilibrium Equations and Boundary Conditions 

Equilibrium equations and boundary conditions that are work conjugate to the strains 
appearing in equations (15) are obtained by applying the principle of virtual work. The statement 
of this principle for the shells considered herein is given by 



& 7 {/\ A,A 2 d^ 1 d§ 2 



5%' E A,A 2 d^ 1 d§2 + 



6%^ ds 


(30) 


where <S70 , is the virtual work of the internal stresses and 6 7i/ ^ is the virtual work of the external 
surface tractions acting at each point of the shell reference surface A depicted in figure 1. The 

B 

symbol &%( represents the virtual work of the external tractions acting on the boundary curve 

d A that encloses the region A, as shown in figure 1. The specific form of the internal virtual work 
needed to obtain the equilibrium equations is obtained by substituting equations (26) and (28) into 
equation (19b). The result of this substitution yields 
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5^ = ({^} T [do] + {^} T [k 0 ]){5u} 



( 31 ) 


The pointwise external virtual work of the tangential surface tractions q, and q 2 and the normal 
surface traction q 3 is 


The surface tractions q p q 2 , and q, are defined to be positive-valued in the positive ^ 2 -, and in- 
coordinate directions, respectively, as shown in figure 2, and include the effects of a live normal- 
pressure field (see Appendix A), denoted by the superscript "L" and dead surface tractions, 
denoted by the superscript "D. " The boundary integral in equation (30) represents the virtual work 
of forces per unit length that are applied to the boundary dA of the region A, and it is implied 
that the integrand is evaluated on the boundary. The symbol ds denotes the boundary differential 
arc-length coordinate, which is traversed in accordance with the surface divergence theorem of 

Calculus. For many practical cases, the domain of the surface A is given by < ^, < b , and 

a 2 < < b 2 , and the boundary curve dA consists of four smooth arcs given by the constant values 

of the coordinates and § 2 , as depicted in figure 2. The general form of the boundary integral 
is given by 



(32a) 


where 



(32b) 
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(33a) 



[5^(N-a,) + 6^(N 



ds 


where 


6‘^'k, = N,Su, + S,6u 2 + Q,6u 3 + M,8q), + M 12 6cp 2 


(33b) 


= S 2 8u, + N 2 5u 2 + Q 2 8u 3 + IVL^cp, + M 2 6cp 2 


(33c) 


In equations (33); N p S p Q, , N 2 , S 2 , and Q 2 are the c; and components of the external forces 
per unit length that are applied normal, tangential, and transverse to the given edge, respectively, 
as shown in figure 2. Likewise, M, and M, are the components of the moment per unit length with 
an axis of rotation that is parallel to the given edge, at the given point of the boundary. In addition, 
M 12 and M 2 j are applied twisting moment per unit length with an axis of rotation that is 
perpendicular to the given edge, at the given point of the boundary. 

The equilibrium equations and boundary conditions are obtained by applying "integration- 
by-parts" formulas, obtained by specialization of the surface divergence theorem, to the first 

integral in equation (30). For two arbitrary differentiable functions f(i„i 2 ) and g(£„i 2 ), the 
integration-by-parts formulas are given in general form by 



(34a) 


(34b) 


where N is the outward unit-magnitude vector field that is perpendicular to points of dA, and 
that lies in the corresponding reference- surface tangent plane. In addition, a ] and a 2 are unit- 
magnitude vector fields that are tangent to the and coordinate curves, respectively, at every 
point of A and 3 A, as shown in figure 1. 

The integration-by-parts formulas are easily extended to a useful vector form by noting that 

T 

the product {v} {w} represent a linear combination of scalar pairs. Thus, the vector forms of 
equations (34) are given by 


33 


f f {f} T ^d|,d§ 2 = - f f ^i|L{g}d^d^+ f ^^(N-2,) ds (35a) 

J J A J J A Jd A 


ii 




f f ^-{g}d§,d§2+ f ^ (N^jds (35b) 

J J A J d A 


Applying these equations to the left-hand side of equation (30), and using equation (31) for b7{/ { 
gives the following results: 


JY '({*>>'] +w , [k,])^ 


d{6u} 


ai. 


ii 


d 

w, 


A 2 ({^}‘[di] + {^} 1 [ k 1 ])){ 8u / d ?i d §2 

+ j ({^} T [<ix] + {^} T [ k 1 ]){ Su } (N • a,) ds 

Jd A 


(36) 




1 d{5u) 

^^A,A2 d ^2 = 


{^} T [d 2 ] + {^} T [k 2 ]) A ^ 

A,({^} T [d 2 ] + {^}'[k 2 ])J{5u}d§ 1 d§ 2 

+ j ({^} T [d 2 ] + {3£} T [k 2 ]){& u } (n • a 2 ) ds 

Jd A 


ii. 


d 


(37) 


Jf 


A.W'lkujAJ-1 


d(6u} 


d §i d i 2 


ii 


d 

1 1 d 


A, d§, 


(A 2 {^f} T [k n ]) {6u} di,d^ 2 


( 38 ) 



__1 

A,A 2 d§, 


(A J {»} , [k 11 ]){5u} - )(n • a.) ds 


34 



35 


JI 


A 2 {?, 


a t3{6y°} 
J 3|, 


d§,d§ 2 = 


I I i( A 2{5 7 i} T ){ 6 Y°} d §!d§ 2 + I {? 1 } T {8 y°} (N • aj ds 

J J A ' JdA 


( 42 ) 


JI 


t3{6y°} 

a '<*> - V 


d=,dS. = 


I | ^^( A i{5 7 2} T ){Sy°} di,d§ 2 + I {? 2 } T {5 y°} (N • a 2 ) ds 

J J A 2 JflA 


(43) 


From these equations, it follows that 


JI 


JI 


b7{/ l A,A 2 d^,d § 2 = { 2 ?} {6u} + { 2 ?} { 6 y 0 } d§,d| 2 + I bTO, ds ( 44 ) 


I 

•/ a a 


where 


2?} T = A.A 2 ({^} T [do] + {^} T [ k o]) 


a 

' 3|; 
d 

d "%2 


a.(<*> t [«i,] + {*}>.]) + J + x“^( A '{^> T )[ k .d 


A, d £; 2 


(45a) 


f 

A >( W T N + {*}>,]) + ^-^(A 2 {%} T )[k 12 ] + ^(A l{ ^} T )[k 22 ] 


A, d % 2 


{2 ?} T = a,a 2 {2} t - ^(a 2 {?,} t ) - ^-(a,{? 2 } t ) 


(45b) 


Next, equations (2) are used to get 
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P22 Mil 


A jA 2 


A, d§, 


[<*>>.]]- 


A 2 


+ {W} [k„] 

(W’M 


1 9 


( 


A,A2^, 
1 9 


A i {*}'[k,] + ^j|^(A I {*> T )[k 11 ] 


+ 


V 


j_J_ 

a, a § 2 


A,{^} T )[k 12 ] 


A,A 2 ai ; 


A i {^}> 2 ] + -f^(A 2 {^} T )[k 12 ] + iA( AlW T )[k 22 ] 


l 






1 / y5l I _/M T . 1 /~r T 1 /c?x T 1 3 /c7\ t 1 d /c? \ n 

* ,27, - ,2, + ,?.) - p J7,, - A|5?1 ^> -A.s?,' 7 -' 


1 a 


A ]A 2 


The boundary integral appearing in equation (44) is given by 


I ds = I [5^(N*a,) + 5^ I2 (N*a 2 ) 

J dA J a a 


ds 


where 


m, = 


{*} [«!,] + {«} [k.] 


{5u} 


AjA; 


a 


A 2 {^}' , [k 11 ]) + -^(A,{^} T [k 12 ]) 


< 6u } 


, . t f n 1 a{ 5 u) , . t f 1 1 a{ 5 u) , .t, „ , 

-<*> [ k 4x;-^7 L -W + (M 


b7t/‘ = ({*> , [d.] + {*> T [k, ]){8u} 


+ 


1 

[ 3 f 

A ,A 2 



a 

a§ 2 




a, *, <*>>•> 


Next, equations (32a) and (44) are substituted into equation (30) to get 


(46a) 


(46b) 


(47) 


(48a) 


(48b) 
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2 ?} T - AjA 2 {q} T ){5u} + (2 ?} T {&Y°} 



[hTV* - «&£) ds = 0 


(49) 


as the statement of the principle of virtual work. Because {8u} and {5y 0 } are independent 

virtual displacements, that are generally nonzero, the localization lemma for the Calculus of 
Variations yields the equilibrium equations 

x L<2?> T -{q> T =!0} (50a) 

{2?} T =(0} (50b) 

In order to identify the strain-displacement relations of several different theories that appear 
often in the literature, equation (16a) is expressed as 


/ 

\ 



and equations (4) as 



e,-, -T 


cp 


(?!, §2) - 2 C = 


r(tPi + c 2 cp 2 

+ 2 C| 

(® 11 ) 

2 

+ eu 

e°2 

+ 2cp] 

r(tP 2 + c 2 cp 2 

+ 2 C| 

( e 22) 

+ e [2 

, e > 2 

-2<p] 

+ cp(Cp 2 + c 3 

e ii(e° 2 - cp 

+ e 22 ( 

0 

e i2 

+ cp) 

„ \ C 3 U t 

1 3u 3 




^1- R, ' 

A.ai. 




c 3 u 2 

1 du 3 





^ r 3 - 

A, d% 2 





1 5U 2 1 

du, 

, U ‘ 

\ 

+ U2 



A.ai. a 

2 

Pi. 

P 22 J 




n A 


(51) 


(52a) 


(52b) 


(52c) 


With this notation, equations ( 1 6a) and (8) are recovered by specifying c, = c,, = c 3 = 1 . Specifying 
c, = 0 and c, = c 3 = 1 gives Sanders’ strain-displacement relations, and specifying c, = 0, c, = 0, 
and c 3 = 1 gives Sanders’ strain-displacement relations with nonlinear rotations about the 
reference-surface normal neglected. In addition, specifying c, = c 2 = c 3 = 0 gives the strain- 
displacement relations of the Donnell-Mushtari- Vlasov theory. Accordingly, 
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{&0 = 



^ (l + c,e°, )5e°! + c 1 (e ° 2 + cp) 6 e ° 2 + cpjbcp, + (c,e ° 2 + c,cpj5(p 

(l + c 1 e 2 2 ) 5 e ° 2 + c,(e ° 2 - qpjSe^ + qp 2 6 cp 2 - ( c^e”, - c 2 cp) 8 cp > 

t (53) 

c,(ei 2 - tp) 6 e I | l + c,(ei 2 + cp) 6 e 22 + ^2 + c^e',) + e 22 ) J 6 e ”2 
+ tpjScpj + cpj6(p 2 + - e^jbcp 


where 




c 3 5u 3 1 <38u 3 

r, 


(54a) 


Scp 2 (^i, £ 2 ) = 


c 3 6u 2 1 36u 3 


r 2 a 2 


(54b) 


6(p(?„ i 2 ) = jC 


i d5u, i d6u, 8u, Su, 

— 1 H - 

A, (3^j A 2 d§ 2 p ] i P 22 


(54c) 


Substituting equations (23a)-(23c) and (54) into (53), the matrices in equations (26) are expressed 
as 


[d»] = 


3 cp, c,(l + c 3 )e° 2 + (c, + c 2 c 3 )cp 


R, 


2p, 


,(l - c 3 )e° 2 + (c 2 c 3 - c,)cp l+ Cl e“ 2 


2p ii p 2 2 

c 3 cp 2 2 + c,( 1 - c 3 )e°, + c,( 1 + c,)e° 2 c,(e° 2 + cp) 


R, 


+ ■ 


2pj 


,(c 3 - l)e° 2 + (c 2 c 3 - c,)cp 1 


+ c,e, 


2P22 Pll 

c 3 qp 2 c,( 1 + c 3 )e° 2 - (c, + c 2 c 3 )tp 

R 2 2p„ 

c 3 tp, 2 + c,( 1 + c 3 )e°, + c,( 1 - c 3 )e 22 c,(e° 2 -cp) 


1 + c,e°, 

R, 

1 + c,e 22 

R 2 

e° 2 - cp e° 2 + cp 


R, 


2p 2: 


R, 


+ 


R, 


(55a) 
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1 + cfi n 

C j / 0 \ C / ° 

■y( e i2 cp j + -y( c i e i2 

+ c 2 cp) 

- CP! 

[d,] = 

0 

y(e° 2 - cp) - y(c,e ° 2 

- c 2 cp) 

0 


Cj(e ° 2 - cp) 

1 + ^ ( 1 _ c 3 )e n + -^( 

+ c 3 )e 22 

- cp 2 


y(e : 2 + cp) 

- y(c,e “ 2 + c 2 qp) 

0 

0 


tK _(P ) 

+ y(c,e “ 2 - c 2 qp) 

-t 0 

1 + c je 22 

-cp 2 


1 + y(l + C, 

)e°, + y(l - c 3 )e 22 

-i(e ° 2 + cp) 

-cp, 


Moreover, 


[«.] - J-[4.] + J-[4.] = 

M 22 Mil 

c 3 <Pi + c '( e ‘- + jp) _ 1 + C|e", 
R-l P 1 1 P 22 

c .( e ° 2 -tp) | i + c, e ; 

P 1 1 P 22 

c 3 cp 2 2 + c ] (e,, + e 22 ) 2c,cp 
R 1 P 1 1 P 22 


_ c i(e‘i’ 2 + cp) _ 1 + c,e 21 
P 22 Pn 

c 3 (p 2 _ c 1 ( e 12 ~ cp) + 1 + c,e 22 

^■2 P 22 P 1 1 

c 3 cp, 2 + c,(e n + e 22 ) 2 c,cp 

^2 P22 P 1 1 


1 + Cje n qpj 

Ri P22 

1 + C t C 2 2 — 2 

1^2 Pll 



e t2 + qp j ( gp 

R 2 p 2 


In addition, substituting equations (54) in to equation (22b) yields the revised matrices 


[k„] = c. 


1 


A.ai.lR, 
1 


P ] |R-2 


P22R, 


1 a 


( 1 A 


A 2 d^; 


V R W 


1 


1 


1 

+ — 


1 a 


2p n 1 R, RJ A j a§, 



+ 


Ri 


0 

0 

0 


(55b) 


(55c) 


(56) 


2 

Pn 


(57a) 
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Cj_ 

R. 


0 


0 0 


0 


C3 ( A 

2\R, 



0 

1 

P22 

1 


(57b) 


[kj = 


R, 


0 


C 3 

2lR, 


1 


1 

Pn 

0 

1 

P22 


(57c) 


that appear in equation (28). Equations (29d)-(29f) remain unaltered. By applying these revised 
matrices to equations (46) and (50), the equilibrium equations are found to be 


1 d3?n + 1 d3?!2 
A, 3 |, A 2 3§ 2 


1 s^ l2 | 9^2 

A, 9 ?=, A 2 9§ 2 


2^,2 ^22 c 3 Q 13 c 3 3 / 1 

Pn P22 R i 2 A 2 9 § 2 [ -{R, 

7 'tw ~ Tin 171 n C 3 Q 23 C 3 3 

^ — rv r - 

P11 P22 ^-2 2 A. j 




+ P, + q, - 0 


(58a) 


J +P 2 + q 2 = 0 (58b) 


1 9Q 13 + 1 9Q 23 Ql 3 Q 23 ^11 ^22 

A t A 2 9^ 2 P22 Pn R 1 R-2 


(58c) 


?2i ?i 1 1 9p n 1 9p 2 i 

PlI P 22 9^i -^2 9^ 2 


(58d) 


?22 “Pl2 1 ri?, 1 9p 2 2 

P11 P22 A, 9 ^, A 2 3§ 2 


(58e) 


where 

O B + _L^Ii + Wu -^22 _ »o 

A, 9 §, A 2 3^2 p 22 Pn 


(59a) 
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(59b) 


~ ^ 1 97X i2 + i d^ 22 + + 2%r, 2 

A, d§, A 2 d% 2 p n p 22 


P,= 


[*1.<P. + 2 A + ^^22 ) J ^ [^11®11 "*" ^12( e i2 *p) 

— — ^ n (e 12 + (p) + ^22(^12 — *p) ”^12 ( e 1 1 + ^22) — ~ [^22^22 — ^n®n + 2^ 12 cpj (59c) 

0 r (^,i -^22)cp + 2(^22e; + ^i2e;,)] 


2 a 2 a^ 


^2 = -^ 


[^22^2 ^12*Pl j 2 A ag [(^ n "*" ^ 22)^1 "*■ . ^22^22 ^ 12(^12 "*" ?*) 


2 A, a§ 


A 2 a§ 2 


+ 


^11(^12 + tp) + ^22(^12 cpj + ^^(c 1 1 + 622) [^22^22 ^n®ii ■!■ 2 ^? 12 (pl ( 59 d) 

*22 L J P n L J 


+ 


2A,a§ 1 L 


— ”^22 j^P + 2('% 11 e 12 + ^12622 


^ 3 — r^n^Pi + ^12^2] [^n*Pi "*■ ^12^2] . -,t- “I” ^22^2] 

A| wsi L P 22 ^ A 2 di 5 2 L 

H |^12*Pl "I" ^22*pbJ ~ j^ll^ll "I" ^ 12(^12 — *p) j — ^22^22 + ^12^12 ) j 


(59e) 


Before the boundary conditions can be obtained, the boundary integral given by equation (43) 
must be reduced further. In particular, by noting that 


it is seen that the integrals 


j: 

j: 


(n • a, )ds = A, d^ 2 

(60a) 

(n • a 2 )ds = - A, d^, 

(60b) 

{»} f )E ,'( N ' a ') ds 

(61a) 

{»} T [ kja -^ M ( N. Sj ) ds 

(61b) 
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can be integrated by parts further, by using the product rule of differentiation, to get 

^} T [k J 

J dA 



1 d{bu\ j - ^ x 

J — hr 1 N* a, 
A 2 d% 2 1 lJ 


a, a|, 1 2j 


ds = 


Using this result, equation (47) is expressed as 


I ft?#* ds = I • Z 2 ) ds + ({^f} T [k 12 ]{8u}) 

J dA J 5 A 


(63) 


where 


^> T [k 12 ]{5u}) 3A - f ^|-[({%} T [k I2 ]){6u}]d§ 2 - f ^|-[({^) T [k 12 ]){5u} 
rt aa rt a a 


d|, (64a) 


= 


*} T [<*.] + {»}>.] + x;i( <*>’I k "]) ( 8u > 


1 

a 

A iA 2 

31. 


-(A 1 {»} T [k„]) + 4-( Al {»} , [k„]) 


{Su} (64b) 




43 


w l2 = 


W T [d 2 ] + {«} T [k a ] + 


(5u) 


+ 


1 

a 

A jA 2 

ag, 


A 2 {*r}‘[k 12 ]) + ^r( A > W'tM) {6u} (64c) 


<*}> + <?=}'{**'} 


On an edge given by = constant; dg, = 0, (n • a, ) = 1 , and (n • a 2 ) = 0. For this case, 


f 6^-ds= r 

J dA J a. 


8^ n A 2 ]g i=COMt!mt d| 2 + f({3fc} [ k 12 ] {^u} ^ 


§ ! = constant 


(65a) 


where 


{VC) [k 12 ]{5u} /Si=constmt 


= [{«} [ k ,z]{Su> 


g ! = constant 


Likewise, equation (33a) has the form 


I b7(/ E ds = I 

J dA J a. 


b 2 


b7V, A:\. . . dg 2 


g 1 = constant 


(65b) 


(66a) 


with 


= 


N,(g 2 ) + 


3M.(g 2 ) 

R, 


Su, + 


s,0M + 


3 M 12 (g 2 ) 

R, 


6u 2 + Q,(g 2 )Su 3 


M 1 (g 2 )36u 3 M 12 (g 2 )a6u 3 


a, agj a 2 ag 2 

where equations (54) have been used for the virtual rotations. Noting that 


(66b) 


r 

, . a6u , 

t r 

[dM p (g,L 

1 

m, 2 (? 2 ) ^ 

dg 2 = - 

g 1 = constant 

•/ a. 

J 6u 3 

dg 2 


g 1 = constant 


dg 2 + (M 12 (^)6u 3Si= _)" 2 (67) 
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equations (66) are expressed as 


I' a^f[ 

dA mJ a. 


m' E ds - d?, - (M.ai.lSu 


3 ^ , = constant 


(68a) 


with 




N,fe) + 


CjM,(^ 2 ) 

R, 


8u, + 


s. 0 M + 


C 3 M 12 (^2 


R, 


8u, 


+ 


Qi(M + 


dM 12 (^ 2 ) 

d§. 


Su 3 - 


M,(g 2 )d6u 3 

Aj d%, 


(68b) 


The matrix form of equation (68b) is given by 


with 


In addition, 


5^ = {^} T {8u} - {^ 1 } T [k 11 ] 


1 d{5u) 

A, d%, 


f 

N,+ 


CjMj 


R, 


A 


f \ 
\*if 


s,+ 


c 3 M 12 

R 


Qi + 

l 



y 



(M 12 (§ 2 )6u 35i = con5tan ,)" : =( 


{m,} [k,,]{5u| 


^ i = constant 


(69) 


(70a) 


(70b) 


(71) 


Thus, enforcing the boundary integral term in equation (49) for an edge given by = constant 
by using equations (64b), (65), and (67) and applying the localization lemma of the Calculus of 
Variations yields 
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m, - w, = 




+ 


A ,A 2 


M 


<6u) 

(A,{W}> 11 ]) + ^-(A l <»} T [kj) 
.}’ - {w> T )[k„] + {?,}'{&/■) = 0 


d 

Wl’ 


I 611 ) 


(72a) 


(» 


> T -{*,> T )[k„]{8»} 


= 0 


§ , = constant , 


(72b) 


1 d{6u} 

where ( 8u } , — — — — , and {fry 0 } are arbitrary virtual displacements. The component form of 

A, d§, 

these equations yield the following boundary conditions for the edge given by ^ = constant: 


%„(l +c,e‘; i )+% 12 c 1 (e‘; 2 -(p) + ^ 11 ^ = N 1 (| 2 )+M 1 (^ 2 )^ or 6u, = 0 (73a) 


Ttn + + ^ 22 )^ -y ^n(2e 12 + (pj — ^? 22 cp + 2^ 12 e 2 


+ ^ 12 ^ 


v r 2 rJ l(?z) 


or 6 u 2 = 0 (73b) 


Ri 


— 1 / x . dM 12 (^ 2 ) 

Qis + "^ — 1 - (^n<Pi + ^12^2) = Qi (^2) H ^ — — or 6 u 3 = 0 


(73c) 


«.i = M,(i 2 ) or 


I dSu 3 
A, d§, 


= 0 


(73d) 


?n=0 or 6 y°3 = 0 


(73e) 


7,2 = 0 or SyL = 0 


(73f) 


and 


2^12 = M 12 (^ 2 ) or §u 3 = 0 


(73g) 


at the corners given by = a 2 and '% 2 = b 2 . Examination of equations (20c) reveals that 7, and 
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5^ appearing in these boundary conditions correspond to forces per unit length associated with 
through-the-thickness distributions of a n and o l2 , respectively, that supress transverse-shearing 
deformations of the plate edge face. 


On an edge given by = constant; d§, = 0, (n • a,) = 0 , and (n • a 2 ) = 1. For this case. 


i 


/»b, 


h7V x ds = 


+ (({*}>„]{««}) 


= constant 


c/ a. 


(74a) 


where 


({^} T [k 12 ]{6u})^ coM ) = ([{^} 1 [k 12 ]{5u}] 


J §2 = constant J 


(74b) 


Likewise, equation (33 a) has the form 


6^> = 

c/ dA •J a. 


b, 


6^ B ds=| 5 ^ e , A, I d|, 

h ■ 1 t2 1 = constant 


(75a) 


with 


m, = 


sMA + 


c 3 M 2 ,(^|) 

R. 


6u, + 


N 2 (S,) + - 


3 m 2 (^) 


R, 


6 u 2 + Q 2 (?,)6u 3 


M 2 (g,)a6u 3 M 21 (g 2 )a6u 3 

A 2 d§2 A, 


(75b) 


where equations (53) have been used for the virtual rotations. Noting that 


1 


b, 


. ,d6u, 

M,, fe. — 4 


di, = - 


1 J §2 = constant 


1 


b, 


dM 21 (^j) 

— ^ 8l L 

d|, 


| 2 = constant 


d|. + (M 2I (^ 1 )6u 31 _ constant )" 1 (76) 


equations (75) are expressed as 


w>= [«*UA,] 

«/ dA i/a. 


§ 2 = constant 


d§, - (m 21 (^,)6u 


3^2 = constant 


(77a) 
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with 


(77b) 


m, = 


s 2 (?,) + 


3 M 21 (1 1 )' 


R, 


6u, + 


n 2 (1,) + 

Q 2 (^.) + 


c 3 M 2 (|,) ~ 

R 2 

dM4|i| 

d£, 


5u, 


8u, 


M 2 (Ii) d6u 3 

A, d§ 2 


The matrix form of equation (77b) is given by 



- {% 2 } T [k 22 ] 


1 d { 6u } 

A 2 d| 2 


with 


c , C3M-21 

2 R, 

lj/ \ = J XT . C 3^2 

\*'-f ~ \ ^ + R 
dM 2 , 


Q 2 + 


d?, 


(78) 


(79a) 




(79b) 


In addition, 


(M 21 (? 1 )8u 3S2=con5tant )" i , = ( 


{Vt 2 } 


[ k , 2 ]{8u} 


^2 = constant 


(80) 


Thus, enforcing the boundary integral term in equation (49) for an edge given by = constant 
by using equations (64c), (74), and (76) and applying the localization lemma of the Calculus of 
Variations yields 


m n - m E2 = 


2 T 

tl) 


< 6U ) 


(*} T [dJ + W> 2 ] + 

-(A 2 {^} T [k 12 ]) + -^(A,{^} T [k 22 ]) {8u} (81a) 


1 

a 

A ,A 2 

3|, 


- {^} T )[k 22 ]^-^^ + {? 2 }>Y°} 
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T )[k 12 ]{8u} 


= 0 


'§2 = constant 


(81b) 


1 a ( 6u } 


where {8u}, — — tz — and {Sy"} are arbitrary virtual displacements. The component form of 

A. 2 ^S2 

these equations yield the following boundary conditions for the edge given by = constant: 


- y(^i. + ^)<P + y (^11 - ^ 22 )<P + 2(v 22 e° r _ + ^e",) 


or Su, = 0 (82a) 


% 22 (l +c 1 e;) + ^ 12 c 1 (e; 2 + ( P )+^ 2 2^ = N 2 (| 1 ) + M 2 (^ 1 )^ 0 r Su 2 = 0 (82b) 


Q 23 + 


1 d%? l2 / % , 

~7 ZZ [^12*Pl ^^22 T 2 ) — 02(^1 

Aj oc, j 


dMJ?,) 

+ ^ - or 8u 3 = 0 (82c) 


and 


, , 1 d6u 3 

^22 = M 2 (?,) or 

-^2 ^S2 

?ii =0 or 8 y°3 = 0 
?, 2 = 0 or by ° 23 = 0 

^i 2 = M 21 (i 1 ) or 8u 3 = 0 


(82d) 

(82e) 

(82f) 

(B2g) 


at the corners given by = a t and = b r Examination of equations (20d) reveals that "P l2 and 
7 2 , appearing in these boundary conditions correspond to forces per unit length associated with 
through-the-thickness distributions of a 12 and o 22 , respectively, that supress transverse-shearing 
deformations of the plate edge face. 


Alternate Form of the Boundary Conditions 

In the present derivation of the boundary conditions, and those given by Koiter, 6,9 ' 10 the 
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derivatives 


^ dSu 3 
A, dlj 


and 


2 36u 3 
A, d^ 2 


have been taken as the basic displacement parameters 


along the edges = constant and t;, = constant, respectively. In contrast, Sanders and Budiansky 5 
8,11 use the virtual rotations hep and hep, as the corresponding basic displacement parameters. For 
these basic displacement parameters, equation (73a) and (73d) are replaced with 


^n^l + Cje,, ) + %i 2 Ci(e 12 cpj — N, (^ 2 ) or bu , — 0 
= M,(i 2 ) or hep, = 0 

for the edge = constant. Similarly, equation (82b) and (82d) are replaced with 


(83a) 

(83b) 


"^ 22 ( 1 ^i^ 2 2 j ^ 12 ^ 1(^12 tp ) — bJ 2 (^,) or bu 2 — 0 (84a) 

^22 = M 2 (§,) or hep, = 0 (84b) 

As pointed out by Koiter 9 (see part 3, p.49), these alternate boundary conditions are completely 
equivalent to those presented herein previously, and are completely acceptable. 


Thermoelastic Constitutive Equations for Elastic Shells 

Up to this point in the present study, the analysis has a very high fidelity within the 
presumptions of "small" strains and "moderate" rotations. The constitutive equations are 
approximate in nature and, as a result, are the best place to introduce approximations in a shell 
theory. The constitutive equations used in the present study are those for a shell made of one or 
more layers of linear elastic, specially orthotropic materials that are in a state of plane stress. 
These equations, referred to the shell ('§,, gj coordinate system are given by 


o n 

o 22 

0,2 


QnQnQ .6 

Ql2 Q22 Q 2 6 

Q ,6 Q 26 Q 66 




(85a) 


/ \ 
\ / 


c 55 c 45 / y,3 \ 

C 45 C 44 \ Y23 / 


(85b) 


The Q y terms are the transformed, reduced (plane stress) stiffnesses of classical laminated- shell 
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and laminated-plate theories, and the C- terms are the stiffnesses of a generally orthotropic solid. 
Both the Q ;j and the Ck terms are generally functions of the through-the-thickness coordinate, 

for a laminated shell. The GL terms are the corresponding transformed coefficients of thermal 
expansion, and 0(§ p | 2 , | 3 ) is a function that describes the pointwise change in temperature from a 
uniform reference state. An in-depth description of these quantities is found in references 226 and 
227. 


The work-conjugate stress resultants defined by equations (20) are the only stress resultants 
that appear in the virtual work, equilibrium equations, and boundary conditions. Thus, the shell 
constitutive equations are obtained by substituting equations (15) into (85), and then substituting 
the result into equations (20). This process gives the two-dimensional shell constitutive equations 
by the general form 


w 


W 


{?.} 

= 



{2} 



[c 00 ] [C 01 ] [C J [Co,] 

[Cj [cj [c 12 ] [c 13 ] 
[c 20 ] [cj [c 22 ] [c 23 ] 
[cj [c 31 ] [c 32 ] [cj 
[cj [cj [C 42 ] [cj 


[C 04 ] 

[Cm] 

[c 24 ] 

[CJ 

[C 44 ] + [c J 


/p°\ 
\ £ / 


/ Y n 

\X / 


1 9 

kW, 

1 9 
A 2 9^2 


</} 

{y} 


/ v °\ 

\Y / 


\0«/ 

If* \ 

\0,/ 

{0 2 } 

If* \ 

\©3/ 

{e 4 } 


( 86 ) 


where 


[C J = 


h \ l + R. 

2 


1 + 


R, 


[S.]' 


Ql, Q,2 Ql6 

Q12 Q22 Q26 

Q.6 Q 2 6 Q« 


[Sj]d§ 3 for i,j£{0, 1, 2, 3, 4} (87a) 


[Css] = 


h \ l + R, 

2 


1 + 


R. 


[Ss] 1 


c c 

v -"55 ^45 

c c 

'—'45 '—'44 


[S s ]d§ 3 


(87b) 
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( 88 ) 


/<a \ _ 
\O k f - 


[SJ 


Q„ Q,:Q„ 
Ql2 Q 22 Q 2 , 
Ql6 Q 26 Q ft 



©(ip §22 § 3 )d^3 


Substituting equations (17) into equation (87) and using the shorthand notation z, = 1 + 
z 2 = 1 + — , and Z = z, + z, + y (z 2 - z, ) yields the following exact expressions 


la. 

R, ’ 


[Co.] 


J 


[C.,] = 




r 





^Qn 

z, 

Q 12 



Q,2 

-Q 22 

Z 2 

Z _ 

2^ Q “ 

h 

2 

z _ 

~ Ql6 

2Zj 

Z _ 

2z)° 26 

z 2 _ 

7 Q« 

4ZjZ 2 

Z 

- Qn 

Q 12 

z, + z 2 

Z 

2z, 


Ql2 

Z| Q 
z 2 ^ 2 

z, + z 2 


2z 2 


di. 


i:* 5 * 


Z\ Z 2 ry ^ 

A ^ ^< 6 < 

4z,z 2 


§ 3 di 3 


(89a) 


(89b) 



_ h 

2 

-F.^jQ,, 

Z, 

— F 2 (§ 3 )Q,6 

[c 02 ] = 


F.fejQ* 

F 2 (| 3 )Q2 6 


h 

2 

Zf,( S ,)Q,. 

^&)Q« 


d|; 5 


(89c) 
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j 


h 

2 



F,(^)Q I6 

F 2 (| 3 )Q,2 


-F,(| 3 )Q 26 

-F 2 (| 3 )Q22 


Z 2 

Z 2 


^F.(l3)Q« 

1^)0. 

/_h 

2 






(Q 12 
\ P 22 

+ 22 Ql6 If, 

z, Pi, j 

LQh, + Z 
\ P 22 Z 


Z| Q22 _j_ Q26 |p 

Z 2 P 22 P 11 


2 \ Z 2 P 22 Z lPl 


d§3 


z, CL + Q A 


z, p 22 


2 v Z 'P, 


+ 


Pn y 

CL 

Z 2 P 22 


f 2 

f 2 


d§ 3 



1 ”* 2 

^Qn 

z, 

Q,2 

z, + z, _ 
2z 


Q,2 

—0 
7 V 22 

^2 

z ‘ + Z: q 26 

2z, V26 


z, + z 2 

Z, + Zt 

O 

(z,+z 2 )^ 


2z Ql6 

2z 2 

4z,z 2 




z 2 


Z, 


F,(^)Q„ — F 2 (^ 3 )Q„ 


F.(^)Q,2 F ; (i 3 )Q:, 


§ 3^3 


Z, + Z, Z| + Zo 

-F,(^)Q 16 -d— F 2 (?,)Q a 


2z, 


2z, 


(89d) 


(89e) 


(89f) 


(89g) 
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2 


' 


F,(i 3 )Qi 6 

F 2 (i 3 )Q,2 



^F 2 (^ 3 )Q 22 


Z 2 

Z 2 

h 

z, + z 2 , 

— F(S. )Q« 

z, + z, , 

2z, ' F *-' Q -| 


Ssdk 


(89h) 



r»- 

2 


Qn + ^oJ 

P22 z i Pn ; 


Z 2 P22 Pn ; 

— f^Q + — Q 

2z,z 2 P22 Pn 


'Q.+SQ^f, 


3 P22 


^Q* + 

Z 2 P22 


Z 1 Pll 

Q ,2 




^ Q “ + ^ Ql 

2 z 3 Z 2 \ P22 Pn 


(89i) 


[C 22 ] 



f, 2 (^)Qh f 1 (§ 3 )f 2 (§ 3 )q 16 ^ 
f,(^)k(? 3 )q I6 f 2 (^ 3 )q 66 z, 3 


(89j) 




F, 2 (^ 3 )Q 16 F.(| 3 )F I (i 3 )Q 12 

f 1 (§ 3 )f 2 (^ 3 )q 66 f;(§ 3 )q 26 


d§ 3 


(89k) 


[C»] = 


J _ h 
^ 2 


Z, p 


On 

P 22 ^ 1 h 1 1 y 

Q26 + Z 2 Q 6 6 Ip p 

P 22 Z 1 P 1 1 / 


( Ql, . Z 2 Q 11 

f,f 2 

\ P22 

z i Pn y 

( Q 6 6 . 22 0,6 1 

f 2 2 

y P22 

z i Pn ) 


d?3 


(891) 
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[C„] = 


F, 2 (^,)Q 66 f 1 (§ 3 )f 2 (§ 3 )Q 2 , 
f,(^ 3 )f 2 (§ 3 )q 26 f;(^)q 22 


j-\ 


z 2 


d§ 3 


(89m) 




[C M ] = 


J _ h 

47 2 


Q66 _j_ Q 26 .p- 

P li Z 2 p 2 2 


Q,, 

Pn 


+ 

z 2 P22 j 


F,F, 


'Q,. + giQ. A 


Pn 

Q12 

Pn 


P22 j 
Z 2 P22 j 


f,f 2 


d§ 3 


(89n) 


fj 


[Cm] = 


./ -h - 


z 2 Q 66 , 2 Q 2 , 


+ 

z i P|j P11P22 z 2 


+ 5^% |F, 2 


2 P22 


Z 3 Ql6 , Ql2 + Q66 . Z,Q, 


Pn 


P 1 1 P 22 


J2 P22 


^2 Q11 

'■ Pi 


Ql2 + Q66 Z, Q 26 jp p 


P 1 1 P 22 


gaQu t 2 Q, 6 

Zip^ P11P22 Z 


J2 P22 J 

+ ^%W 

2 P22 / 


d§ 3 (89o) 




[C»] = 


•M 


— ( z,F/ — A I C, 


z,F/ - 


R, 


z 2 F 2 ' 


z,F/ 


R 


z ^'-f;P 


Ac 

R 2 14 


Z 2 


Z|, z 2 F 2 '- a] c 4 


r 2 


d§ 3 


(89p) 


and [CJ = [Cij] for i, j E {0, 1, 2, 3, 4}. Next, the thermal stress resultants { © k ) are obtained 
by first expressing the temperature-change field by 


i 3 ) = 0(l.. i 2 )G(i 3 ) 


(90) 


such that 


{e k } = e(i 1 ,i 2 ){ 0 t } 


(91) 


with 
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r 2 




Qn Qi: Qi6 

Ql 2 Q22 Q26 

q 16 q 26 q 66 


I- 

s 


a 

«22 }G(^ 3 )d§ 3 


kG {0,1, 2, 3, 4} (92) 


Specific expressions for {® k } are obtained by substituting equations (17) into equation (92). 

The constitutive equations given by equations (86), (89), and (92) include the effects of 
transverse-shearing deformations in a very general manner. Specifically, the constitutive 
equations of a given transverse-shear-deformation theory depend on the choices for the functions 
F ] (§ 3 ) and F 2 0§ 3 ), which are required to satisfy F^O) = F 2 (0) = 0 and F,'(0) = F 2 '(0) = 1. With these 
functions specified, expressions for the transverse shearing stresses are then obtained by 
substituting equations (4e) and (4f) into equation (85b) to obtain 


fVOM 1 + 


R, 


F.fe) 

R, ^ 


F 2 '(0 1 + 


1 + 


G55Y 13 + 


R, 


R, 


1 + 


fife) 

R ^ o 

C 45 y 2 


(93a) 


R 3 


F/fe) 1 + 


Oot 


R, 


F,(i 


^45 Y 13 ^ fc 


FM 


(93b) 


1 + 


R, 


1 + 


R, 


In refined transverse-shear-deformation theories, it is desireable, but not always necessary, to 
specify F^,) and F,(| 3 ) such that a 13 = o 23 = 0 at = ± h/2. For example, two similar choices for 
the pair F,(^ 3 ) and F 2 (^ 3 ) will yield different stress predictions yet their stiffnesses, obtained from 
equations (89), will yield nearly identical predictions of overall buckling and vibration responses. 
In contrast, for accurate stress analyses, one might expect Fj(§ 3 ) and F 2 (^ 3 ) to account for the 
inhomogeneity found in a general laminated-composite wall construction. Examples of various 
choices for F 3 (§ 3 ) and F 2 (^ 3 ) that have been used in the analysis of plates and shells are found in 
references 193-195. Some examples are discussed subsequently. 

Constitutive equations that correspond to a first-order transverse-shear-deformation theory 
are obtained by specifying 

F ,(^ 3 ) = F 2 (? 3 )=?3 ( 94 a ) 


For this fundamental, first-approximation case 
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+ 


(94b) 


cx, 3 = 


C55Y 13 C 45 y 2 


1 + 


R, 


1 + 


R, 


C 45 y 13 

H 


1 + 


R, 


1 + 

R, 


(94c) 


Thus, the functions specified by equation (94a) lead to expressions for the transverse-shearing 
stresses that do not satisfy the traction-free boundary conditions o 13 = cr 23 = 0 at q, = ± h/2. 
Moreover, these choices for the distibution of the transverse-shearing stresses and strains, given 
by equation (94a), do not reflect the inhomogenous through-the-thickness nature of laminated- 
composite and sandwich shells. 

A more robust transverse-shear-deformation theory is obtained by specifying 

F,(l 3 ) = F a (l 3 ) = ^3 - (95a) 


Upon substituting these functions into equations (3) and neglecting the nonlinearities, it is seen 
that the displacements U, and U 2 are cubic functions of q, and that U 3 has no dependence on 
at all. As a result of this character, a transverse-shear-deformation theory based on equation (95a) 
is referred to herein as a {3, 0} shear-deformation theory. For this case 


CToq 


1 - 

1 + 2 
3rJ 

(*)■ 


1 + ^ 



R, 




/ _ x 2 

1 - 


m 


L 3R.J 

l h ) 


C 5SY 13 + 


1 + 


3R, 


1 + 


21, 

h / _ 

C 45 y 2 


1 - 


1 + 


I*. 

R, 


c 45 y I3 + 


1 + 


R, 


3 R 2 


2 |j 

h 


1 + 


C 44 y 2 


R. 


(95b) 


(95c) 


The functions specified by equation (95a) also lead to expressions for the transverse-shearing 
stresses that do not satisfy the traction-free boundary conditions a 13 = cr 23 = 0 at q, = ± h/2 unless 
h/R and h/R, are negligible. If 6/Z? denotes that maximum shell thickness divided by the 
minimum principal radius of curvature, then 


/°»\ 

\°23 /?,= * 





(° \ 

c c 

^55 '-'45 

< 

Y 13 \ 

c c 

^45 44 

?3=±| 

l Yz3 / 


(96) 
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For moderately thick shells with <£/;? = 0.1, the coefficient in equation (96) involving /l/IR is equal 
to 0.035. This result suggests that the error in the traction-free boundary conditions may be within 
the error of the constitutive equations for the practical range of values 61 <0.1. However, these 
choices for the distibution of the transverse-shearing stresses and strains, given by equation (95a), 
also do not reflect the inhomogenous nature of laminated-composite and sandwich shells. 

Expressions for Ffe.) and F,fe) that satisfy o 13 = a 23 = 0 at = ± h/2 are obtained by 
assuming cubic polynomials for F,(| 3 ) and F 2 fe), with a total of eight unknown constants, and 
then enforcing F,(0) = F,(0) = 0 and F,'(0) = F 2 '(0) = 1 in addition to the four traction boundary 
conditions. This approach yields 





3 1 h 


1- JJJl 
12 R. 




(97a) 


Fjy- 


1 + 1 

If) 

few 

(f) 

l-lkl 

(if 




(97b) 


However, R, and R, are generally functions of ('§,, | 2 ), which violates the requirement that F, = 
F.fe) and F 2 = F,('§.). If the terms involving h/R, and h/R, are neglected, then equation (95a) is 
obtained. A simple choice for F^) and F 2 fe) that satisfies a 13 = a 23 = 0 at f = ± h/2 and avoids 
the presence of R, and R 2 is given by 


Fids) = F 2 (§ 3 ) = | 3 - + j4&) 5 (98) 

h h 

These functions and their first derivatives vanish at = ± h/2. 

A highly refined transverse-shear-deformation theory for laminated-composite and sandwich 
shells is obtained by specifying Ffe :i ) and F 2 (|,) based on "zigzag" kinematics. This approach 
has been presented in reference 228 for laminated-composite and sandwich plates. This approach 
satisfies the traction-free boundary conditions a 13 = a 23 = 0 at t;, = ± h/2, and yields functional 
forms for F,(^ 3 ) and F 2 (§ 3 ) that account through-the-thickness inhomogeneities. Moreover, for a 
general inhomogeneous shell wall, F,(| 3 ) and F 2 (| 3 ) are found to be different functions. 


Simplified Constitutive Equations for Elastic Shells 

The exact forms of equations (87) and (88) are obtained by integrating equations (89) and 
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(92) exactly, once the through-the-thickness distribution of the elements of the constitutive matrix 
in equation (85) are known. Generally, this process leads to very complicated functional 
expressions for the stiffnesses defined by equations (89). However, for laminated-composite 
materials, the elements of the constitutive matrix in equation (85) are modeled as piecewise- 
constant functions and the integration of equations (87) and (88) poses no problems. 

In the present study, the approximate nature of the constitutive equations is exploited to 

P P 

simplify the shell constitutive equations by expanding the functions of ' and ’ appearing in 

K, K 2 

P 

equations (89) in power series and then neglecting terms that are third-order and higher in — and 

p 

’ . In addition, it is noted that 

Ko 


- 1 + 


z 2 


( 1 

1 \ 

- 1 . x 2 

(ri 

rJ 

- ^(Ss) 


» 1 + 


{ \ l\ 

, 1 , , 2 

( r ^ _ rt ) 

^3 - ^(§ 3 ) 


i-i + i -L--L 
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Following this process yields the following constitutive equations for a laminated-composite 
shell: 
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In these equations, x = 1 gives the constitutive equations that are second order in — and 3 

R, R 2 

that are generalizations of those put forth by Fliigge 171 for isotropic shells. Setting x = 0 gives 
constitutive equations that are first order. It is also important to note that the matrices obtained 
from equation (101) for k = 0, 1, and 2 correspond to the [A], [B], and [D] matrices, respectively, 
of classical laminated-plate theory (see reference 226). Similarly, the thermal parts of the 
constitutive equations are given by 


W = 



+ 


hu 

R 2 

1 

h^ 

R, 


K 


2(R, + R 2 ) 


1 

R 


1 

R, 


0 

0 

2 

h [2 


(e.) = 


+ 


h~, 

r 2 

2 

h 22 

R. 

2 

h 12 


2(R, +R 2 ) 


(107a) 


(107b) 


(0 \ _ / gn\ + 1 /gn\ 


{0 ) _ / g.T\ 1 /g!'\ 


(107c) 


(107d) 


66 


(107e) 


( e <)=7r 


10 g , 2 

g„ + — 

r, 


20 

gl 2 


gn 

R, 



where 


h„ 

hi 

k 

h„ 




Qn Q,2 Q„ 
Q12 Q22 Q:, 

Q , 6 Q 26 Q a 


«22 ^(§3) G(i-,)dS 
[a, ' 


(108) 





Qn Q ,2 Q„ 

Q,2 Q22 Q„ 
Qi 6 Qn, Q,„ 


laj 


F j (? 3 )(? 3 ) G(i 3 )dS 


(109) 


By examining equations (94a) and (95a), it is seen that the maximum magnitudes of 
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and 


: h 3 are typically greater than zero and less than unity. Using this information, the constitutive 

equations given by equations (93) are simplified, using binomial expansions of the denomenators, 
to obtain 
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Special Cases of the Constitutive Equations 

In the original shell theories of Sanders, Budiansky, and Koiter; the constitutive equations 
used are the simplified, first-approximation constitutive equations of classical Love-Kirchhoff 

linear shell theory. Their constitutive equations follow from neglecting ^ and in 
expressions for the effective stress resultants given by equations (20) when they are used with 
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equations (85) to determine the constitutive equations. A similar set of constitutive equations are 
obtained from equations (106) by neglecting all terms involving principal radii of curvature R 
and R,. In addition, terms involving the radii of geodesic curvature p u and p 22 are neglected. These 
terms originally entered the constitutive equations through the matrix [SJ in equations (87). To 
see the rationale for neglecting these terms, this matrix is expressed as 
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Next, the fact that the maximum magnitudes of ‘ h ' and h are typically greater than zero 

and less than unity is used again. In addition, the radii of geodesic curvature measure the bending 
of the coordinate curves within the tangent plane at a given point of the reference surface. 

Typically, p u and p 22 are substantially larger than the shell thickness h and, as a result, h and 

Pu 

h have very small relative magnitudes. Therefore, it follows that constitutive matrices based on 

P22 

[SJ can be neglected based on the inherent error in the constitutive equations. Based on the same 
reasoning, equations (110) are approximated as 
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A detailed derivation of these equations is presented in Appendix B. Additionally, constitutive 
equations are presented in Appendix B for transverse-shear deformation theories that include a 
first-order theory, a {3, 0} theory, and a zigzag theory. 


Effects of "Small" Initial Geometric Imperfections 

The effects of "small" initial geometric imperfections are obtained in the present study by 
following the approach presented in reference 221. With regard to the strain fields, the effects of 
"small" initial geometric imperfections appear only in the nonlinear membrane strains given by 
equations (5). These effects are obtained by replacing the normal displacement u 3 (i,, '§,) with 
u 3 (§„ |,) + w'(!„ '§,), where w'(| p | 2 ) is a known, measured or assumed, distribution of reference- 
surface deviations along a vector normal to the reference surface at a given point. Then, all terms 
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involving w'd,, | 2 ) appearing in a given membrane strain, that correspond to an unloaded state, 
are subtracted from that given strain. Applying this process to equation (51) yields 
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Equation (26) is then expressed as 
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where [d 0 ], [dj, and [d 2 ] are given by equations (55) and where 
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As a result of the "small" initial geometric imperfections, three additional terms appear in the 
counterpart of equation (46a); these terms are 
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Including these additional terms in equations (46) and (50) yields the equilibrium equations 
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The terms in equation (122d) arise from the live pressure load. 

Additional terms associated with the initial geometric imperfections also appear in the 

T 

boundary conditions given by equations (72) and (81). In particular, {3?} [d,] in equation (72a) 
is replaced with {3?} [d, + d,] , and [d 2 ] in equation (81a) is replaced with 
{7t) [d 2 + d 2 j. As a result, equations (73a)-(73c) become 
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Likewise, equations (82a)-(82c) become 
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In addition, the alternate boundary conditions given by equation (83a) becomes 
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and that given by equation (84a) becomes 
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Expressions for the displacements U,, U 2 , and U 3 are obtained by replacing u 3 d„ 1,) with 
u 3 d„ 1 2 ) + w'd,, 1 2 ) in equations (52) and (3), and then eliminating terms involving w‘d,, 1 2 ) in 
equation (3) that are left over when u,, u 2 , and u 3 are set equal to zero. This process gives 
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A complete resume of these fundamental equations is given in Appendix C. 
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Concluding Remarks 


A detailed exposition on a refined nonlinear shell theory that is suitable for nonlinear limit- 
point buckling analyses of practical laminated-composite aerospace structures has been 
presented. This shell theory includes the classical nonlinear shell theory attributed to Leonard, 
Sanders, Koiter, and Budiansky as an explicit proper subset that is obtained directly by neglecting 
all quantities associated with higher-order effects such as transverse-shearing deformation. This 
approach has been used in order to leverage the exisiting experience base and to make the theory 
attractive to industry. The formalism of general tensors has been avoided in order to expose the 
details needed to fully understand and use the theory in a process leading ultimately to vehicle 
certification. 

The shell theory presented is constructed around a set of strain-displacement relations that are 
based on "small" strains and "moderate" rotations. No shell-thinness approximations involving 
the ratio of the maximum thickness to the minimum radius of curvature were used and, as a result, 
the strain-displacement relations are exact within the presumptions of "small" strains and 
"moderate" rotations. To faciliate physical insight, these strain-displacement relations have been 
presented in terms of the linear reference- surface strains, rotations, and changes in curvature and 
twist that appear in the classical "best" first- approximation linear shell theory attributed to 
Sanders, Koiter, and Budiansky. The effects of transverse-shearing deformations are included in 
the strain-displacement relations and kinematic equations, in a very general manner, by using 
analyst-defined functions to describe the through-the-thickness distributions of transverse- 
shearing strains. This approach yields a wide range of flexibility to the analyst when confronted 
with new structural configurations and the need to analyze both global and local response 
phenomena, and it enables a consistent building-block approach to analysis. The three- 
dimensional elasticity form of the internal virtual work has been used to obtain the symmetrical 
effective stress resulants that appear in classical nonlinear shell theory attributed to Leonard, 
Sanders, Koiter, and Budiansky. The principle of virtual work, including "live" pressure effects, 
and the surface divergence theorem were used to obtain the nonlinear equilibrium equations and 
boundary conditions. 

A general set of thermoelastic constitutive equations for laminated-composite shells have 
been derived without using any shell-thinness approximations. Acknowledging the approximate 
nature of constitutive equations, simplified forms and special cases that may be useful in practice 
have also been discussed. These special cases span a hierarchy of accuracy that ranges from that 
of first-order transverse-shear deformation theory, to that of a shear-deformation theory with 
parabolic through-the-thickness distributions for the transverse-shearing stresses, and to that 
which includes the use of layerwise zigzag kinematics without introducing additional unknown 
response functions into the formulation of the boundary-value problem. In addition, the effects of 
shell-thinness approximations on the constitutive equations have been presented. It is noteworthy 
that none of the shell-thinness approximations appear outside of the constitutive equations. 
Furthermore, the effects of "small" initial geometric imperfections have been introduced in a 
relatively simple manner to obtain a nonlinear shell theory suitable for studying the nonlinear 
limit-point response. The equations of this theory include tracers that are useful in assessing many 
approximations that appear in the literature. For convenience, a resume of the fundamental 
equations of the theory are given in an appendix. Overall, a hierarchy of nonlinear shell theories 
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have been presented in a detailed and unified manner that is amenable to the prediction of global 
and local responses and to the development of generic design technology. 
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Figure 1 . Coordinate system and unit-magnitude base-vector fields for points 
of undeformed shell. 
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Appendix A - Live Normal Pressure Loads 


For a live normal pressure load p , the pressure p depends on the deformation of the shell 
reference surface. Thus, 


p = q 3 (?i + u„ + u 2 , + u 3 ) 

Expanding p in Taylor Series and retaining terms up to first order gives 
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The differential force due to the live pressure acting on a deformed- shell reference surface is given 
by 

dF = pd§£ (A3) 


where «(ii> i 2 ) is the unit-magnitude vector field normal to the deformed reference surface and 

d§ is the deformed image of the reference- surface differential area dS. The surface area dS is 
given in principal-curvature coordinates as 

dS = A,A 2 d^d^ 2 (A4) 

For "small" strains and "moderate" rotations, the analysis of reference 221 indicates that 

d§£ = ds(qp,a, + cp,a 2 + ( 1 + e? t + e 22 )nj (A5) 


where a,(§,,^ 2 ), a 2 (§ l5 § 2 ), and n(|„ | 2 ) are the unit-magnitude base- vector fields associated 
with points of the undeformed reference surface, cp, and cp, are the rotations defined by equations 
(8), and e°, and e 22 are the linear deformation parameters defined by equations (7). Substituting 
equations (A2) and (A5) into equation (A3) gives 
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where terms involving products of displacements, strains, or rotations are presumed negligible. 
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Appendix B - Special Cases of the Constitutive Equations 


Several special cases of the constitutive equations are presented subsequently in which all 
terms involving the principal radii of curvature are neglected. Likewise, terms involving the radii 
of geodesic curvature are also neglected. This approach leads to equations that are consistent with 
the classical shell theory and the shell theories of Leonard, 7 Sanders, 8 Koiter, 910 and Budiansky. 11 
Thus, equations (85) reduce to 
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In addition, equations (92) become 
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With these simplifications, equation (86) reduces to 
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In addition, the constitutive equations for the transverse-shearing stresses are approximated as 
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(B22d) 
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It is important to note that for Fj(§ 3 ) = F 2 (§ 3 ), equations (B22) reduce to 
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where ? l2 = ? 2l . 

First-Order Shear-Deformation Theory 

The constitutive equations of a first-order shear-deformation theory are obtained by setting 
F i ( 1 3 ) = F-,(g 3 ) = J§ 3 in equations (87)-(94). In particular, equations (87) and (88) become 
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where the right-hand-sides are given by equation (86). Likewise, equation (91) becomes 
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Moreover, equation (94) becomes 
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where the right-hand- side is given by equation (93). Applying these simplifications to equations 
(B23) yields 
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As a result of 'p l2 = “P 2V Substitution of equations (B27) and (B28) into equations (B23) reveals 
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that e p n = 7fi n , 7 22 = and ? 12 = Thus, the constitutive equations are expressed in terms of 
the nomenclature of classical shell theory as 
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where k 44 , k 45 , and k 55 are transverse- shear correction factors that are used to compensate for the 
fact that the transverse-shearing stresses are uniformly distributed across the shell thickness and, 
hence, do not vanish on the bounding shell surfaces. Additionally, the two equilibrium equations 
given by (58d) and (58e) become 
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Applying these matrices to equation (20e) and using equations (13) gives 


R>\ = 
\ 2 ™) 



(B37) 


For the special case in which the stiffnesses and the thermal coefficients appearing in 
equations (85) are symmetric through the thickness, the following additional simplifications to the 
constitutive equations are obtained 
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{3, 0} Shear-Deformation Theory 

The constitutive equations of a {3, 0} shear-deformation theory are obtained by setting 

F 1 (i 3 ) = F 2 (i s ) = ^ 3 - A (^) 3 (B40) 


in equations (87)-(91) and (94). These functions yield parabolic distributions of transverse- 
shearing stresses across the shell thickness. Substituting these functions into equation (B22d) 
reveals that the transverse- shearing stresses vanish at = ± h/2. Next, using equation (B40), 
equations (96) and (97) become 
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where the right-hand-sides are given by equation (95). Likewise, equation (100) becomes 
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and equation (103) becomes 
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where the right-hand-side is given by equation (102). Applying these simplifications to equations 
(B23) gives 
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For the special case in which the stiffnesess and the thermal coefficients appearing in 
equations (85) are symmetric through the thickness, the following additional simplifications to the 
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constitutive equations are obtained 
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Shear-Deformation Theory Based on Zigzag Kinematics 

When the transverse-shearing stresses are approximated by equation (B22d), the analysis 
presented in reference 228 for laminated-composite plates is directly applicable to the present 
study. In particular, the functions F,(^,) and F 2 (| 3 ) are expressed as 

F 1 (i,) = f 1 (i J ) + vr , (i I ) (B50a) 

F 2 (£,) = f 2 (i,)+'I r f , (S.) (B50b) 

where the ply number n G { 1 , 2, N} and N is the total number of plies in the laminated wall. 
The functions f,(| 3 ) and f,(?= 3 ) are continuous functions with continuous derivatives that are 
required to satisfy 
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(B51c) 


f /(*|) = f 2 '( ± |)=° (B51d) 

for all values of -h/2 < < h/2. The functions Tf 1 '(!,) and T, ’(§,) are referred to in reference 

228 as "zigzag enrichment functions" that account for shell wall inhomogeneity and asymmetry. 
These functions are given by 
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In equations (B52) and (B54), the superscript M denotes the ply that occupies the shell reference 
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surface, c;, = 0. The derivatives of F,(| 3 ) and F 2 (§ 3 ) are 
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For a homogeneous shell wall, the zig-zag enrichment terms in equations (B50) and (B56) vanish. 
Thus, the derivatives of the functions f,(^ 3 ) and f,(| 3 ) represent distribution of transverse-shearing 
stresses in a homogeneous shell wall. For the parabolic distribution of transverse shearing stresses 
commonly found in the technical literature for homogeneous plates, 
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The shell-wall stiffnesses and thermal coefficients appearing in equations (B22) are obtained 
by substituting equations (B50), (B56), and (B58) into equations (101)-(109) and performing the 
through-the-thickness integrations. The resulting expressions are lengthy and are not presented 
herein. 
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Appendix C - Resume of the Fundamental Equations 

The fundamental equations needed to characterize the nonlinear behavior of shells with 
"small" initial geometric imperfections are presented in this appendix. The parameters c v c 2 , and 
c 3 appear in these equations and are used herein to identify other well-known shell theories that 
are contained within the equations of the present study as special cases. In particular, specifying 
F,(§ 3 ) = F,0§ 3 ) = 0, neglecting the initial geometric imperfections, and setting c t =c 2 = c, = 1 gives 
the nonlinear shell theories of Budiansky" and Koiter . 9,10 Similarly, specifying c 1 = 0 and c 2 = c 3 
= 1 gives Sanders’ nonlinear shell theory, and specifying c t = 0, c, = 0, and c 3 = 1 gives Sanders’ 
nonlinear shell theory with nonlinear rotations about the reference-surface normal neglected. 
Furthermore, specifying c t = c 2 = c 3 = 0 gives the Donnell-Mushtari-Vlasov 192 nonlinear shell 
theory. 

Displacements and Strain-Displacement Relations 

The fundamental unknown fields in the present study are the reference-surface tangential 
displacements Ujd,, | 2 ) and u,(|,, | 2 ), the normal displacement u 3 (|„ and the transverse- 

shearing strains and Y 23 (^ 1 * ^)- The corresponding displacements of a material point 

(Ip l 2 > 13) are § iven b y 
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where w'(i,, | 2 ) is a known field that describes the initial geometric imperfections in the 

unloaded state. The functions F,(^,) and F,(t.) are user-defined and specifiy the through-the- 
thickness distributions of the transverse-shear strains. These two functions are required to satisfy 
F^O) = F 2 ( 0 ) = 0 and F/( 0 ) = F,'( 0 ) = 1. The nonzero strains at any point of the shell are given by 
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where the reference-surface membrane strains are given by 
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+ C|e,,— op,- — 

R, A,^, 


(C9) 


e 22 (i,, ? 2 ) = e 22 + ^-(cp 2 + C 2 cp 2 ) +ic, (e 22 ) + e° 2 (e° 2 - 2cp) 


0 w ] 1 dw' 

+ C 1 Cp?— — — 

R 2 T A 2 a^ 2 


(CIO) 


Y° 2 (?i- S 2 ) = 2e, 2 + cp,cp 2 + c e;(e; - cp) + e 22 (e° 2 + cp) 


c,w 


qW, „ Y C,w , a , 1 dw‘ 1 3w‘ 

+ T 7 ie ' 2 “ "’) + ir( e ” + <*>) - 'F'xiair ‘ 1 : 


(Cll) 


with the linear deformation parameters 
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(Cl 2) 


1 3u i _ + }h 

A,^, p„ R, 




<4(1, . i,) = J- + i- + £ 

A 2 ^^2 P 22 R -2 


2 ^ 12^19 §2) “ 


1 3Uj 1 du 2 Uj u 2 


A 2 5§ 2 Ajd^i p n p 2 


XnMO = 


1 dcp, cp 2 


p„ 


xlMO = 


1 dcp 2 qpj 


A 2 ^^2 P2 


^2) _ 


1 dqp 2 , cp i aqp! cp 2 1 1 


+ + 


- <P 


A| Pn A, d£ 2 P 22 l Ri R 


and the linear rotation parameters 


/ fc fc \ _ c 3 u i 1 au 3 

(Pl (^ 1 ’^ 2 ) R| A,a|, 


<p 2 (ii. £ 2 ) = 


c,u 2 1 3u 3 


r 2 A 2 a| 2 


<p 


s 2 ) - 


1 ( 1 du 2 1 <3u, u, u 


+ — - + 


“ \ A j <3§! A 2 d^ 2 Pn P 22 


Equilibrium Equations and Boundary Conditions 

The equilibrium equations are given by 


l d% u | d7l n 2 7t n 7t u - 7t 22 c 3 Q 13 

+ — — + + 

1 1 


A, d§, A 2 d§ 2 p, 

+ 


c 3 d 


2A 2 d§ 2 


R, R, 


+ + q t + + q, - 0 


(C13) 
(Cl 4) 
(C15) 
(Cl 6) 
(Cl 7) 

(Cl 8) 
(Cl 9) 
(C20) 


(C21) 
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(C22) 


j_d7t_ 

A, d§, A 2 d'E , 2 


n + 1 ^ 


Tty 


+ 


2^ 

+ 


12 + C aQ23 


+ 


Pll 

c 3 


R, 


2A[ 


R 


1 

R 


+ + q, + ^ + q 2 — 0 


1 dQ n i dQ 23 Q 13 Q 23 Ti u Ti 22 i i 

A, A 2 dl 2 p 22 p n R, R 2 43 4 


9n _ 9u _ J_d9n _ 1 d? 21 
Pn P22 A, A 2 dE , 2 

“Pll "P \2 1 d “?, 2 1 3^22 

Pii P22 A, A 2 3^ 2 


where 

q 1 &7x n t 1 t m n -m 12 2M n 
A, d§! A 2 d ^ 2 P22 P11 


Q 


23 


1 9^12 1 93*22 -^22 2^,2 

Aj A 2 Pn P22 



L 1 <3w‘ 


L 1 dw' 

q, A 2 d% 2 


\ 

V 

/ 


L 

q 3 


w w 

R, Ri 


A 


dq 3 i 

+ — ~~w 

d§ 3 


J 


Py = 


[*ii<Pi + ^12^2] 2 ')t [ ^ 1 1 ^22) ■*■ ^ ^ j^ii^n + ^12(^12 q) 

- ~^n( e i2 + q) + % 22 (e° 12 - cp) + 7t u (t u + e 22 )| - y~[% 2 c° 22 ~ + 2^ 12 q 

d 

Cl 1 (7t n -^22)^ + 2(^22612 + ^ 12 e°,) 


2 A 2 d§ 2 


(C23) 

(C24) 

(C25) 

(C26) 

(C27) 

(C28) 

(C29) 
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|^ 22 (p 2 ^12*Pl j 2 A [(^ n ”^22 ) ^P j "*" . -,g j%22&22 "*" ^1 2 (®12 + T ) 


2A,a§, 


A,a§ 2 


+ ^11(^12 + Cp) + ^22(^12 tpj + ^12(611 + 622) [^22^22 ^ll^ll + 2 '% 12 cpj (C 30 ) 

P22 L J p n L J 

(^11 — ^22)*? + 2(^11612 + ^12^22)! 


2Aj a^.L 


P 3 = - 


1 3 

A, 31 , 


(*,l«Pl "*” ^12*?2 j [«llCPi “I” 2 J . ->s- ^^l 2 Tl + ^?22*P2 j 


a 2 3§ 2 


l^^^CPi + ‘%22*p 2 j ^11^11 + ^(eis *p)j 1^22^22 ^ 12(^12 + 


<P 


(C31) 


Pl = ir 


1 3w‘ 1 3w' 

^11”! ^ *" ^12 


A, 3£, ,2 A 2 3§ 

Cj 6 


c, 3 


A 2 a§: 


W 

% % 


A,a§, 

c,w‘ 


« w_ 
?11 R, 

■% u ^22 

R, R 2 


c,w' / 1 1 

r^ + r; 


(C32) 




1 dw' 1 dw 1 

^^19 + ^£99 

A, 31, A, 3§ 2 


C[ 5 


A,3§, 


w 

* a R, 


c, 3 

n, w ' 

c,w' 

f«n 

% 2 \ 

A 2 3§ 2 

^ 22 r 2 

Pn 

[r, 

R 2 j 




(C33) 




1 


A,3|, 

1 


1 3w‘ 1 dw‘ 

^'a^^aTsi; 


a 2 3|, 


1 3w‘ 1 3w‘ 

* ,! ATai: + *” Al 3|7 


1 3w' 1 dw‘ 

^11 21 1 " ^19 

A,3|, 'A 2 


1 

Pn 


1 3w‘ 1 3w‘ 

^19 1" ^99 “ 

‘A, 3§ t “ A, 3^2 


- c 


% 2 ^ + % L ^_ 

R, R 2 R| R] 


(C34) 


The boundary conditions for an edge given by = constant are given by 

+ 7 g l 2 c,{e ° l2 - qp) = N,(£ 2 ) or u,=D,(i 2 ) 


% 


1 1 u W 

1 + c, e„ + — 


R, 


(C35) 
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C? ( _ _ _ . \ C j 


^12 + + ^22)^ + 


+ ^,2y 


2 

/ 3 1 


3? n (2e° 2 + cp) - ^ 22 cp + 2 1t x 
M I2 (§ 2 ) 


' w 1 ' 


C99 + 


R, 


or u, = D 2 (^ 2 ) 


rI-r; =s,(y+ r 


Ql 3 + 


1 dWn 

a 2 a§ 2 


_ ( 1 dw' 1 ! ( 1 aw 

* 1 , i‘ Pl kW,) + ' 2 V 2 


- Qi(l2j + _ diT" or 


u 3 = D 3 (§ 2 ) 


^r n =M,(^ 2 ) or cp, = <J>,(t; 2 ) 
?,, =0 or Yi 3 = r 1 (§ 2 ) 
7,2 = 0 or Y 2 3 = r 2 (§ 2 ) 


and 


^i 2 = M I2 (1 2 ) or u 3 = D 3 (§ 2 ) 


at the corners given by t;, = a 2 and = b 2 . The boundary conditions for an edge given by 
§2 = constant are given by 


%2 - y(^ll + ^22) 


cp+ 2 

f 3 1 




+ '% 22 (2e" 2 - cp) + 271 x 
M 21 (S.) 


w 


or u, = D , (§ j ) 


R, rJ S ^') + R[ 


71 ,. 


w 


1 + c 1 1 e 22 + ^ 


+ ^i 2 Ci(e° 2 + cp) = N,(|,) or u 2 = D 2 (§,) 


Q 23 + 


l rm i2 


A, a£, 


. , 1 aw‘ \ „ / 1 aw 1 

^12! Ti 7 a - ^ 22 T2 


a 2 as§ 2 


- Q 2 (?,) + ^ |r - or 


u 3 = d 3 (i, ) 


(C 36 ) 


(C 37 ) 

(C 38 ) 

(C 39 ) 

(C 40 ) 

(C 41 ) 

(C 42 ) 

(C 43 ) 

(C 44 ) 
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^22 = M 2 (§ 2 ) or qp 2 = < l > 2(5i) =0 


(C45) 


?2.=0 

or Yi3 = T 1 (§ 1 ) 

(C46) 

?2Z = 0 

or y° 3 = r 2 (§, ) 

(C47) 

and 



W 12 = M 21 (i,) 

or u 3 = D,(§j) 

(C48) 

at the corners given by = a l and = b,. 




Stresses and Constitutive Equations 

The stresses at any point of the shell are given by 



Qll Q,2 Ql6 
Ql2 Q 22 Q 2 6 

Q,a Q 26 Qee 



(C49) 


/ a »\ 

\°®/ 


C55 C 45 / Y 13 \ 

C 45 C44 \ Y 23 / 


(C50) 


where 


©(I.,!*. i s ) = 0 (i..is)G(i,) 


(C51) 


is the known temperature field for material points of the shell. The general constitutive equations 
for a shell are given by 
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i dy° 3 ] 
A, 3?, \ 
i dyh I 
a, 3§, j 



i dy^ 
A 2 a? 2 \ 
i dyh f 
A, dl; 2 J 





JM 

W 




+ 



i d Y ; 3 \ 
A, 3§, \ 
1 dY ° 3 [ 
A, 31;, j 



1 3y° 3 ) 
A, 3? 2 \ 
1 3 yL j 
A, 3? 2 j 


+ 




where the matrices [C.J appearing in these equations are given by equations (89) and the vectors 
{©.i} are given by equation (92). Moreover, [Cj = [C„] . Specialized forms of [C y ] and 
{ ©j } are given in Appendix B. 
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